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Cosmological parameter constraints from SDSS luminous 
red galaxies: a new treatment of large-scale clustering. 
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ABSTRACT 

We apply a new model for the spherically averaged correlation function at large pair 
separations to the meas urement of the clustering of luminous red galaxies (LRGs) 
made from the SDSS bv lCabre fc Gaztanagal (|2009a( ). Our model takes into account 
the form of the BAO peak and the large scale shape of the correlation function. We 
perform a Monte Carlo Markov chain analysis for different combinations of datasets 
and for different parameter sets. When used in combination with a compilation of the 
latest CMB measurements, the LRG clustering and the latest supernovae results give 
constraints on cosmological parameters which are comparable and in remarkably good 
agreement, resolving the tension reported in some studies. The best fitting model in 
the context of a flat, A-CDM cosmology is specified by — 0.261 ± 0.013, fJb = 
0.044 ±0.001, = 0.96 ±0.01, Hq = 71.6 ± 1.2 kms^^ Mpc~^ and erg = 0.80 ± 0.02. 
If we allow the time-independent dark energy equation of state parameter to vary, 
we find results consistent with a cosmological constant at the 5% level using all data 
sets: wde = —0.97 ± 0.05. The large scale structure measurements by themselves 
can constrain the dark energy equation of state parameter to wde = ~1-O5!to i5, 
independently of CMB or supernovae data. We do not find convincing evidence for an 
evolving equation of state. We provide a set of "extended distance priors" that contain 
the most relevant information from the CMB power spectrum and the shape of the 
LRG correlation function which can be used to constrain dark energy models and 
spatial curvature. Our model should provide an accurate description of the clustering 
even in much larger, forthcoming surveys, such as those planned with NASA's JDEM 
or ESA's Euclid mission. 

Key words: cosmological parameters, large scale structure of the universe 



1 INTRODUCTION 

The acoustic peaks imprinted on the temperature power 
spectrum of the cosmic microwave background (CMB) 
have now been measured with impressiv e precision by 
a number of experiments ( Lee et al. I I2001I: iBennet et al. I 



tween some param eters which canno t be broken by CMB 
data alone (e.g. ( Bond et all Il997l : lEfstathiou fc Bond I 
'l999!: lBridle et al. i booC Perhaps the two most important 



20031: iHinshaw et al. II2003I . 



20071 . 12OO9I : I Jones et al. II2OO6I : 



Reichardt et al. II2009D . These observations place tight con- 



straints on the values of many of the fundamental cosmo- 
logical parameters. With the fifth year of integration from 
the WMAP satellite there is seemingly litt le room left for 
any deviation from the bas ic ACDM model l|Dunklev et al. I 
l2009l : lKomatsu et al. However, degeneracies exist be- 
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examples are the curvature of the Universe and the equation 
of state of the dark energy, iude = ftiE/pDE, where Pde is 
the pressure of the dark energy and pde is its density. Mean- 
ingful constraints cannot be obtained on these parameters 
using CMB data in isolation. 

The full potential of the CMB measurements is real- 
ized when these data are combined with other observa- 
tions, such as the Hubble diagram of type la supernovae 
and the larg e scale struc t ure o f the Universe as traced 
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20051 : ISanchez et al. II2OO6I: lAstier et al. II2OO6I: ISeliak et al. I 



2006"; 'Wang & Mukhcricc ' 2006; 'Wood-Vascv ct a l. I I2OO7F 
Spereel et al. 2007: Koma,tsu ct al. _2009j Okum ura et al. 



20081: iXia et al. 



20081 : iFerramacho et al. II2OO9I ) . These com- 



plementary data sets come from the late Universe compared 
with the CMB data, and the interpretation of the observa- 
tions is more complicated and controversial. 

Type la supernovae (SNe) have been proposed as stan- 
dard candles which can probe the luminosity distance - red- 
shift relation. The first strong evidence in support of a cos- 
mological constant ca me from combinin g the SNo data with 
CMB measurements (|Riess et al. I Il998l : jPerlmutter et al. I 
The type la Hubble diagram has come under in- 
tense scrutiny to uncover any hint of non-standardness aris- 
ing from the nature of the host galaxy, po ssible evolution 
with redshift or variation in dust extinction llSuUivan et al. 



20031: [Gallagher et al. II2005I : lEUis et al. II2OO8I : iHowell et al 



20091 ). A recent joint analysis of SNe from different datasets 



suggests that the systematic error on the equation of state 
parameter from a joint CMB and SNo analysis is comparable 
to the size of the random error ( Kowalski et al. . .2008.). 

The power spectrum of galaxy clustering has also been 
used in combination with CMB data (e.g. Percival et al. 
2002). According to the standard lore, the galaxy power 
spectrum on large scales is expected to have a simple re- 
lation to the underlying dark matter spectrum. Moreover, 
the shape of the spectrum is believed to closely follow 
that expected in linear perturbation theory for the mat- 
ter, which can be readily computed given a set of values 
for the cosmological parameters. However, with the avail- 
ability of improving measurements and more refined mod- 
elling of the galaxy power spectrum it has become clear that 
this simple picture is no longer sufficiently accurate to de- 
scribe the data. Additional levels of modelling of the devi- 
ations from linear theory have to be incorporated into the 
analysis. These include empirical models of the nonl inear 
distortion of the power spectrum (jSmith et al. 1 12003 ) and 
possible scale depe ndent bi ases bet ween the clustering of 
galaxies and mass (ICole et al. . ,2005,: iHamann et"ani2008l : 
ICresswell fc Percival Il2009i r Recent studies have cast doubt 
on the accuracy of th ese prescriptions jSanchez fc Cole"! 
I2OO8I : iReid et all 120081 ). In principle, if the simple models 
described the form of the observed power spectrum, then 
the power spectra measured from different galaxy samples 
should yield e quivalent constraints on cosmological parame- 
ters. However. ISanchez fc Colel l|2008l ) found a fundamental 
difference in the shapes of the galaxy power spectra mea- 
sured from the two-degree galaxy redshift survey and the 
Sloan Digital Sky Survey, after attempting to correct the 
measured spectra. These authors found that the red selec- 
tion of the SDSS galaxies results in a strong scale dependent 
bias (see also^Swanson ct al. 2008). Similar scale dependent 
effects have been seen in the power spectrum o f dark matter 
haloes and galaxies mode lled in simulations l|Smith et al. I 
120071 : lAngulo et al. II2OO8I ). 

Recent analyses have not attempted to model the over- 
all shape or amplitude of the power spectrum, as a conse- 
quence of the diSiculties described above in interpreting the 
results from different samples. Instead, attention has shifted 
to a pattern of oscillatory features called the baryonic acous- 
tic oscillations (BAO), which are imprinted on the matter 
power spectrum. These features have been advocated as a 



standard ruler which can be used to measure the distance- 
redshift relation, and hence constrain the dark energ y equa- 
tion of state (iBlake fc Glazebrook . .2003.: ,Hu fc Haiman 



Linder et al I l2003l: ISeo fc Eisenstein I l2003l: IWang 



20od ; lGuzik et al. ll2007l : ISeo fc Eisenstein ll2007l : ISeo et al 
20081 ). The BAO arise from oscillations in the baryon 



photon fluid prior to matter-radiation decoupling. This phe- 
nomenon gives rise to the peaks seen in the power spec- 
trum of temperature fluctuations in the CMB. In the mat- 
ter power spectrum, the oscillations have a much smaller 
amplitude as baryons only account for around 20 per cent 
of the total matter density of the Universe. Furthermore, 
the oscillations in the matter spectrum are out of phase 
with t hose i n the CMB (f S ugivarna iri995l : lEisenstein fc Hu I 



1 19981 . 1 19991 : iMeiksin et al. Ill999l ). The oscillation scale is 



related to the size of the sound horizon at recombination, 
which can be measure d with high accuracy from the CMB 
l|Komatsu et al. Il2009l ). The apparent size of the BAO ruler 
depends upon the parameters woe and Qk, as these deter- 
mine the angular diameter distance out to a given redshift. 
In practice, the BAO are not precisely a standard ruler at 
the level of precision demanded for their interpretation in 
future surveys. Nevertheless, by modelling the appearance 
of the BAO accurately, they are still a valuable probe of 
cosmological parameters (|Sanchez et al. [[20081 : [Smith et al. I 
I2OO81) . Careful simulation work and modelling has shown 
that techniques can be developed which can overcome long- 
wavelength gradients in the power spectrum to yield robust 
constraints on the BAO using the galaxy power spectrum 
(iPerciyal et al. I l2007bl: ISmith et al. I l2007l: IXngulo et al. I 



20081: ICrocce fc Scoccimarro I I2OO8I : iTakahasi et al. I I2OO8I : 
Seo et al. ll20oi v 



The BAO signal has been detected in bot h the 2dF- 
GRS and SDSS surveys (Cole et al. 2005; Eise nstein et al. I 
l2005h . The most powerful BAO measurements currently 
come from sampl es of lu m inous red galaxies (LRGs) 
llCole et al.1 l2005l: iHiitsi I I2OO6I: iPadmanabhan et al. I 



20071: iPer cival et al. '2007^^: 'Okumu ra et al I 120081 : 
Cabre fc Gaztafia ^a 2009a. : Gaztafia ga et al. I l2008bl : 
20091 ). Detections at a lower significance 



Martinez et al. 



have also been reported using galaxy clusters (|Hutsi II2007I : 
lEstrada et al. 1 12009| ). The imprint of these features has 
even been found in th e three point function of LRGs 
((Gaztanaga et al"] l2008al ). Despite this rapid progress, the 
conclusions drawn from measurements o f the B AO remain 
unclear. For example, I Percival et al. I (|2007d ) analysed 
the BAO signal in a joint galaxy sample drawn from the 
Sloan Digital Sky Survey (SDSS) data release five (DR5) 
and the two-degree Field Galaxy Redshift Survey (2dF- 
GRS). Ignoring the information from the amplitude and 
long-wavelength shape of the power spectrum, and simply 
isolating the BAO, the results of Percival et al. showed a 
2.4(T discrepancy with the distance measure ments inferred 
from the supernovae type la (SN) data by lAstier et al. I 
l|2006t ). signaling a possible problem in the modelling of the 
BAO data, or a challenge to the ACDM model. 

Much theoretical work has been devoted to uncover- 
ing scale dependent effects in the BAO an d in improving 



the modelling of the signal in galaxy surveys (^Angul o et al. 
2005; Huff ct al. 2007; Smith ct al. 2007; Angul o et al. 
200&; Smith et al. 12008: 
ISeo et al. ii2008. : .Desiacques 



2008 



Crocce fc Scoccimarro I 
2008h . In order to realize the 
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full potential of the BAO technique as a cosmological probe 
it is essential to quantify any systematics in the signal 
and to understand how the measurements relate to cos- 
mological parameters. First, let us debunk some possible 
preconceptions about BAO. As we have remarked above, 
the BAO are not precisely a standard ruler. In the corre- 
lation function, the Fourier transform of the power spec- 
trum, th e BAO appear as a broad peak at l arge pair sep- 
arations (|Matsubara Il200"3 'l. I Sanchez" et al. I l|2008l ) showed 
that even in linear perturbation theory, the maximum of the 
BAO peak in correlation function does not coincide with the 
sound horizon scale at the percent level of accuracy required 
to fully exploit the measurements expected from forth com- 
ing galaxy surveys. Fu r therm ore, ISmith et al. I (|2008l ) and 
ICrocce fc Scoccimarro I ()2008l ) have shown that both large 
volume numerical simulations and theoretical predictions 
based on renormalized perturbation theory (RPT) indicate 
that the BAO peak in the correlation function is shifted and 
distorted in a non-trivial manner relative to the prediction of 
linear perturbation theory. If unaccounted for, these shifts 
bias the constraints obtained by using the BAO measure- 
ments as a standard ruler. 

Careful modelling of the correlation function is therefore 
required to extract t he cosmological i nform ation encoded in 
large scale structure. ISanchez et al. I (j2008) argued that the 
correlation function is less affected by scale dependent ef- 
fects than the power spectrum and that a simple model for 
the co rrelation function p roposed by Crocc e fc Scoccimarro I 
(|2008l ). based on RPT l|Crocce fc Scocci marro 2006a ij, 
gives an essentially unbiased measurement of the dark en- 
ergy equation of state. This means that information from 
the large scale shape of the correlation function, in addi- 
tion to the form of the BAO peak, can be used to provide 
robust constraints on cosmological parameters. The correla- 
tion function therefore provides a better constraint on the 
distance scale than the more conservative, "BAO only" ap- 
proach required when using the power spectrum (i.e. which 
requires the long wavelength shape information to be dis- 
carded along with the amplitude). 

In this paper we apply this new model to the 
shape of the redshift space correlation function, ^(s), 
of the SPSS Data r elease 6 (DR6) LRGs measured by 
ICabre fc Gaztafiaga I ll2009al'). T he LRG sample analysed 
by Cabre fc Gaztafiaga I ( 2009al l is twice the size of the 
sample used to make the first detection of the BAO by 
lEisenstein et al.l (|2005l ). and has also yielded the first mea- 
surement of the ra dial BAO signal, which constrains the 
Hubble parameter (|Gaztafiaga et aiTlbOOSbl ). We combine 
the LRG clustering information with the latest measure- 
ments of CMB and SNe data. We focus on the constraints 
on the dark energy equation of state and the curvature of 
the Universe, which are the parameters where the extra in- 
formation from the shape of ^(s) can dramatically improve 
upon the CMB only constraints. We also pay special atten- 
tion to the consistency of the results obtained with different 
dataset combinations. 

The outline of the paper is as follows. In Section [2J 
we describe the data used in our parameter estimation. In 
Section [3] we describe the details of our modelling of the 
shape of the redshift space correlation function and com- 
pare it with measurements in N-body simulations. We also 
set out the different parameter spaces that we study and 



describe our methodology for parameter estimation. In Sec- 
tion [3] we assess the impact of the details of the parameter 
estimation technique on the obtained constraints on cosmo- 
logical parameters. Section [5] presents our main results for 
the parameter constraints obtained by comparing theoreti- 
cal models to the CMB data and the correlation function of 
the SDSS-DR6. In Section |6] we focus on the constraints on 
certain distance combinations obtained from the shape of 
^(s). We summarize our conclusions in Section [7] Appendix 
A gives the theoretical motivation for the model we use to 
describe the correlation function and Appendix B gives the 
covariance matrix for the distance constraints. 



2 THE DATASETS 

Here we describe the different datasets that we use to con- 
strain cosmological parameters. The modelling of the cor- 
relation function of luminous red galaxies is described later 
on in Section [3.21 The datasets described below are used in 
different combinations to check the consistency of the con- 
straints returned. 



2.1 The redshift space correlation function of 
SDSS-DR6 LRGs: the monopole 

Luminous red galaxies (LRGs) are an efficient tracer of the 
large scale structure of the Universe. These galaxies are se- 
lected by color and magnitude cuts designed to identify in- 
trinsically red, bright galaxies using SDSS photometry (see 
Eisenstein et al 2001 for a complete description of the color 
cuts). LRGs can be seen out to higher redshifts than galax- 
ies in a simple magnitude limited catalogue, and so map 
a larger volume of the Universe. LRGs have a low space 
density compared to L* galaxies, which means that fewer 
redshifts have to be measured to map out the same volume. 
The low space density, which translates into a higher shot 
noise, is compensated for by the stronger than average clus- 
tering of LRGs, which maintains the signal-to-noise of the 
correlation function at a level which can be measured. 

Here we use the measurement of the 2-point cor- 
relation function of LRGs in redshift space made by 
ICabre fc Gaztafiaga I (|200 9a) shown by the black points in 
Fig. [1] These authors studied the clustering of LRGs in 
Data Release 6 (DR6) of the Sloan Digital Sky Survey 
(SDSS), which has 75 000 LRG galaxies spanning a volume 
of 1 Gpc^ over the redshift interval 0.1 5 < z < 0.47. The 
compa rison between DR6 and the result oflEisenstein et al. I 
|20ok) is plotted in Fig.B4 of ICabre fc Gaztafiaga I (12009^ 1 
and Fig. [T] the two results are in good agreement with the 
DR6 result showing an improve ment of a factor of about 
\f2 in the size of the errors. C abre fc Gaztafiaga I (|2009ah 
carried out an extensive series of tests of their measurement 
of the LRG correlation function. They found that system- 
atic uncertainties in the estimation of the radial selection 
have a small impact on the results (ie see their figure B2). 
The adoption of different weighting schemes also produces 
negligibly small changes (their Fig.B5 and Bll). These care- 
ful tests indicate that this new estimation of the correlation 
function is robust and that any possible systematic effects 
are small compared to the error bars. 

Here we recap two aspects of the analysis of 



© 0000 RAS, MNRAS 000, 000-000 



4 A.G. Sanchez et al. 




U I I I I I I I I I I I I I I I I I I I I I I 

40 60 80 100 120 140 

s(Mpc/h) 

Figure 1. The spherically averaged correlation function of LRGs 
in redshift space. Circles with errorbars show the correlation func- 
tion used in this paper. The shaded region shows the dispersion 
in the correlation function obtained when using the new MAN- 
GLE mask of Swanson et al. (2008) with different completeness 
fractions. The dotted line shows the result for the north stripe 
of DR6. We have also calculated the correlation function for the 
new DR7 (solid line), estimated using a random catalog generated 
from a smoothed version of the selection function. The dashed line 
shows the estimate from DR7 without smoothing. The measure- 
ment from Eisenstein et al. (2005) is shown using triangles; note 
these authors had fewer LRGs and used broader bins. 

ICabre &: Gaztafiaga I ()2009al ) which are particularly relevant 
to this paper: the treatment of the survey mask and the esti- 
mation of the covariance matrix of the correlation function. 

An accurate knowledge of the angular and radial se- 
lection of a galaxy catalogue, including the redshift com- 
pleteness as a function of magnitude and position on the 
sky, is an essential prerequisite for a measurement of clus- 
tering. This information allows the mean density of galaxies 
in the survey to be e stimated. After the release of DR6, 
I Swanson et al. I (|2008l ) provided this information in a read- 
ily usable form, translating the original mask files extracted 
from the NYU Value-Added Galaxy Catalog (Bla nton ct ^ a L. 
20051) , from MANGLE into Healpi x format (jGorski et al. 
19991 ). ICabre fc Gaztanagal |20093) describe how they con- 
structed a survey "mask" for LRGs and tested the impact 
of the mask on clustering measurements using mock cata- 
logues. Using the same techniques, we have also examined 
the correlation function of LRGs in DR7, which has be- 
come available since the submission o f ICabre fc Gaztafiaga I 
l|2009al ). In Fig. [1] we plot a summary of the possible sys- 
tematics in the estimation of the correlation function. The 
estimates corresponding to different versions of the survey 
meisk and different ways of generating a catalogue of random 
points are in remarkably good agreement with one another. 
Furthermore, the results from DR6 and DR7 are in excel- 
lent agreement; DR7 repre sents only a mo dest improvement 
over DR6. The results of iEisenstein et al. I (|2005i ) are plotted 



in Fig. [T] using triangles. Their esti mate is consistent with 
that o f ICabre fc Gaztafiaga ' (2009a') within the errors (note 
that their binning is different, and also, perhaps, their nor- 
malisation, due to a small shift which could be attributed 
to systematics (see Section [2.1.ip ). Whilst there is a mod- 
est improvement in the parameter constraints on using the 
DR7 measurement, in line with the incremental change in 
the number of LRGs and the solid angle covered, we have 
decided to retain the DR6 measurement in this paper, be- 
cause the survey mask has been tested more extensively in 
that case. The DR6 measurement represents a factor of ~ 2 
more volume than that covered by the LRG sample analyzed 

by Eisenstein et al. (2005 ). 

ICabre fc Gaztafiaga I (|2009al ) constructed a covariance 
matrix for the LRG correlation function using mock cata- 
logues drawn from the MareNostrum In stitut de Ciencicis de l 
Espacio (MICEfl) N -body simulations (jFosalba et al. II2008I : 
ICrocce et al. l[2009h . The extremely large volume of this run 
(box size 7680 Mpc) allowed 216 essentially indepen- 
dent LRG DR6 mocks to be extracted. ICabre fc Gaztafiaga I 
(2009 a) investigated different methods to estimate the error 
on the correlation function (see their Appendix A). They 
found that the jackknife error, an internal estimate made 
from the dataset itself (JK; see Norberg et al. 2008), gave 
a reasonable match to the diagonal elements of the covari- 
ance matrix obtained directly from the 216 mock catalogues. 
However, the JK estimate of the off-diagonal elements of 
the covariance matrix is noisier than that obtained from 
the mocks. In this analysis we construct the full covariance 
matrix of the measurement from the correlation matrix es- 
timated from mock catalogues of dark matter haloes with 
similar clustering and abundance to the LRGs, rescaled by 
the JK estimate of the variance from the data, which has the 
advantage of being independent of the cosmological model 
adopted in the N-body simulation. 

2.1.1 The implications of a possible constant systematic 
shift in clustering amplitude 

Small systematic effects including, for example, the integral 
constraint, calibration errors or evolutionary effects, will, if 
unaccounted for, appear as an additive term in the mea- 
sured correlation function: ^(s) — ^truc(s) + ^sys(s), where 
^sys(s) stands for the systematic error. The simplest model 
for Csys(s) is a constant shift, ^sys(s) = K, which we label 
here as the K-shift. Systematic shifts include both effects 
which are unaccounted for in the estimate of the radial selec- 
tion function and angular calibration errors which can move 
galaxies in or out of the galaxy sample. These can introduce 
spurious fluctuations in the observed density. A 1% varia- 
tion in (r — i) color or a 3% shift in r magnitude can intro- 
duce a 10% modulation in the LR G target number density 
lEisenstein et al. Il200il ). However. ISogg et al.1 (|2005l ) found 
that the large scale density variations in the final LRG sam- 
ple are completely consistent with the predictions of biased 
ACDM models, showing that when averaged over large an- 
gular scales, the fluctuations in the photometric calibration 
of the SDSS do not affect significantly the uniformity of the 
LRG sample. Even a spurious number density fluctuation as 

^ http://www.ice.cat/mice 
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large as 5 ~ 5% can only produce a shift K = 5^ < 0.0025. 
As illustrated in Fig[T]tlie potential systematics that we are 
able to check seem to produce shifts of \K\ < 0.001. 

These are small changes, but could our analysis below 
be affected by such a systematic effect, if present? We have 
checked this explicitly by repeating our analysis allowing 
for a constant additive term K as an extra free parame- 
ter, which we marginalize over. We imposed a prior on this 
K-shift as large as \K\ = 0.01. This corresponds to 10% 
density fluctuations, well above the expected systematics in 
the LRG sample. We find that even with this wide prior in 
K, the marginalization over the additive K-shift does not 
change the obtained errors of the cosmological parameters 
and only changes the mean values by less than 0.5a in the 
most extreme cases. This small difference does not justify 
the inclusion of this extra parameter. But note that if we 
want to address the question of what is the absolute good- 
ness of fit, then it might be important to consider such a 
K-shift as an additional degree of freedom. We will come 
back to this question in Sec. [7] 

2.2 The redshift space correlation function of 
SDSS-DR6 LRGs: the radial BAO peak 

The main piece of clustering information we shall use in this 
paper is the monopole or spherical average of the two-point 
redshift space correlation function described in the previous 
subsection. However, there is also useful cosmological infor- 
mation in the form of the correlation function split into bins 
of pair separation p arallel (tt) and perpendicu lar (a) to the 
line of sight, f (a. tt). ^Gaztahaga et al.l (| 2008b^ found a sig- 
nificant detection of a peak along the line-of-sight direction, 
the position and shape of which is consistent with it being 
the baryonic acoustic peak. 

ICaztanaga et al~l (|2008lJ ) measured the position of the 
radial BAO peak (rBAO) using the full LRG DR6 sample, 
and two sub-samples at low, z = 0.15 — 0.30, and high red- 
shifts, z = 0.40 — 0.47 (see their Table II). In this paper 
we make use of these two last measurements, which can be 
treated as independent due to the large separation between 
these redshift intervals. When combined with a measure- 
ment of the sound horizon scale from the CMB, the position 
of the BAO feature can be used as a standard ruler to obtain 
a constraint on the value of the Hubble paramete r, H{z), as 
a function of redshift (Gaztana ga et ahlbOOSbl V Here, in- 
stead of calibrating the radial BAO dist ance rBAo(z) with 
the CM B, we follow the same approach as lGaztafiaga et al. I 
and use the measurement of the dimensionless red- 
shift scale 

AZBAOW = r-BAoW^^^. (1) 

c 

In this equation Hudiz) is the Hubble constant at the 
redshift o f the measurement in th e fiducial cosmology as- 
sumed bv lGaztanaga et al. I l|2008bl ) to convert the observed 
galaxy redshifts into comoving distances. Note that the fac- 
tor H{id{z)/c in Eq. ([ij ensures that Azbao does not depend 
on the fiducial cosmology assumed to obtain the rBAO mea- 
surement. Section [3.31 describes the details of our modelling 
of these measurements. 

Because of the narrow range of a values (0.5 < a < 
5.5 Mpc) used in the radial measurement, the rBAO re- 
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Figure 2. Normalized covariance matrix estimated from 216 
mocks LRG samples. Here we test the covariance between the 
monopole correlation (horizontal axis) and the radial (LOS) cor- 
relation (vertical axis). The correlation is only important on small 
scales and is negligible (less than 10%) for the scales of interest 
here (> 40 Mpc). 

suits are essentially independent from the monopole of the 
two-point function; fewer than 1% of the bins considered in 
the spherically averaged correlation function are used in the 
estimate of the radial correlation. For this reason we treat 
these two datasets as independent. W e tested this explicitly 
using the mock LRG catalogues from lCabre fc GaztafjagaH 
(2009a) (the same set used to estimate the covarince matrix 
in the monopole), to evaluate the covariance of the radial 
and the monopole correlations. Fig. [5] shows the normalized 
covariance: 

„ _ 1 ^ [Uir-Lji)] Mj^-tm 

""M^ -m(i) a.(i) 

where 5m(«)* and ^r(j)'° are the monopole and radial cor- 
relation in the i-th and j-th bins respectivelly measured in 
the fc-th mock catalogue (fc = 1, ...M) and ^ and a are the 
corresponding mean and rms fiuctuations over the M re- 
alizations. As shown in Fig[2] the covariance on the scales 
of interest in our analysis (larger than 40 Mpc) is quite 
small (less than the noise in the estimation, of about 10%) 
which shows that the two datasets are indeed independent 
in practice. This implies that the rBAO measurements ob- 
tained from the radial correlation functions are independent 
from the monopole ^(s). 

2.3 CMB data: temperature and polarization 
power spectra 

The accuracy of recent observations of the CMB mean 
that this is the single most powerful dataset for constrain- 
ing the values of cosmological parameters. The compila- 
tion of CMB measurements we use includes the tempera- 
ture power spectrum in the range 2 < I < 1000 and the 
temperature-polarization power spectrum for 2 < I < 450 
of th e first five years of observations of the W MAP satel- 
lite l|Hinshaw et al. I l2009l : iNolta et al.1 12009| ): the band- 
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power temperature spectrum from the 2008 observations 
of the Arcminute Cosmology Bolometer Array Receiver 
llKuo et al. ] l2007l : iReichardt et al. ] I2OO9I . ACBAR,) over 
the spherical harmonic range 910 < I < 1850; tempera- 
ture and polarization data from the Cosmic Background 
Imager (jReadhead et al. I l2004l . CBI,) with 855 < f < 
1700; observations from the 2003 flight of BOOME RANG 
llRuhl et al. II2003I: iJones et al. |[200alMontrov"eral. ..200& : 
iPiacentiuT et al. I I2OO6I ) in the range 925 < £ < 1400 
and the recent results for the temperature an d polariza- 
tion power spectra measurements from QUaD l|Ade et al. I 
|2008| ) over the range 893 < I < 1864. These measurements 
of the power spectrum of temperature fluctuations in the 
CMB cov er the spherical harmo nic range 2 < I < 1800. 
Following iDunklev et aTl l|2009l ). in order to avoid cross- 
correlations with the WMAP data, we use only the band- 
powers of the small scale CMB experiments that do not 
overlap with the signal-dominated WMAP data. Note that 
the QUaD measurements were not available at the time of 
the "WMAP-f CMB" analysis carried out by Komatsu et al. 
(2008). 



2.4 The Hubble diagram of type la supernovae 

We also consider the constraints provided by the Hubble 
diagram of type la supernovae (SN) as p rovided by the 
UNION sample from ([Kowalski et al. IboOSl ). This compila- 
tion is drawn from 13 indepen dent datasets proc essed using 
the SALT Mght curve fitter (IGuv et al. 1 120051 ) and anal- 
ysed in a uniform way. The sample contains a set of 57 
low-redshift S N, the recent samp les from the SuperNova 
Legacy Survey (|Astier et al. I2OO6I , SNLS,) and the Equation 
of State SupErNovae trace Cosmic Expansion (ESSENCE, 
iMiknaitis et al. 1 12007| ) , t he high-redshift sample from the 
Hubble Space Telescope jR iess et a l. 1|2004 120071) . as well 
older datasets (iRiess et al. 1,1998: jPerlmutter et al. 1 1 19991 : 
iTonrv et al. II2003I : iBarris et al. II2004I ). The final sample is 
the largest available to date, comprising 307 SN which pass 

the selection criteria . 

iKowalski et al. I (|2008l ) suggest a way to include the ef- 
fect of systematic errors when fitting th e SN data. This es- 
timation and those from other authors ( Astier et al. 11200^ : 
IWood-Vasev et al. I l2007l : iHicken et al. II2OO9I ). predict dif- 
ferent systematic errors in the estimated constraints on the 
dark energy equation of state, with values ranging from 5% 
to more than 10%. These approaches differ in the choice of 
which potential sources of systematic errors are taken into 
account and the estimation of their likely magnitude. It is 
our understanding that the community has not reached a 
consensus about the correct way to estimate the effect of the 
systematic errors in the analysis of the SN data. For thi s rea- 
son we follow other authors (e.g. iKomatsu et al. II2009I ) and 
do not include the systematic errors in our constraints on 
cosmological parameters. However, the inclusion of the sys- 
tematic errors in the SN data has important implications for 
the derived values of the dark energy equation of state (see 
Section I5.2|l . This should be borne in mind when compar- 
ing constraints obtained using SN data with those obtained 
using other datasets. 



3 METHODOLOGY 

In this section we summarize the approach used to obtain 
constraints on cosmological parameters. We start by describ- 
ing he different parameter sets that we consider. The para- 
metric model we use to describe the shape of the correla- 
tion function in redshift space is presented 13.21 The the- 
oretical motivation for this parametric form can be found 
in Appendix \K\ In Section [3.21 we also compare the model 
with measurements made using numerical simulations. Sec- 
tion 13.31 describes the model we implement to describe the 
radial BAO measurements. The methodology we follow to 
explore and constrain the parameter spaces is discussed in 
Section 13.41 



3.1 The parameter space 

In this paper we make the basic assumption that the primor- 
dial density fluctuations were adiabatic, Gaussian and had 
a power-law spectrum of Fourier amplitudes, with a negli- 
gible contribution from tensor mod es. From the analys is of 
the fifth year of WMAP data, ^ Komatsu et al. I (|2009| ) did 
not detect any deviation from these hypotheses at the 99% 
confidence limit (CL). Within this framework, a cosmolog- 
ical model can be defined by specifying the values of the 
following eight parameters: 



(fifc, (.Jdm, t^b, T, Us, As, 9, wide). 



(3) 



We now go through the parameters in the above list, defin- 
ing each one and also explaining how the values of other 
parameters are obtained, which we will refer to as derived 
parameters. 

The homogeneous background cosmology is described 
through the various contributions to the mass-energy den- 
sity. These are, in units of the critical density: 0,^, which 
describes the curvature of the universe; ujdm = ^dmh^, the 
density of the dark matter (assumed cold, where h is Hub- 
ble's constant in units of 100 kms~^Mpc~^) and cjb = ^hh^, 
the baryon density. We assume that massive neutrinos make 
no contribution to the mass budget. For most of this paper, 
we assume that the dark energy component has a constant 
equation of state independent of redshift, with the ratio of 
pressure to density given by ii;de. In section [5.31 we relax 
this assumption and analyse models allowing for a time vari- 
ation in this parameter. In this case we use the standard lin- 
ear parametrizatio n given by (jChevallier fc Polarski II2OOII : 
iLinder et aTlbOOSl ) 



woEia) = Wo + Wa{l — a), 



(4) 



where a is the expansion factor and wo and Wa are parame- 
ters. 

The form of the initial fluctuations is described by two 
quantities; the scalar spectral index, Us and the primordial 
amplitude of the scalar fluctuations A^. These parameter 
values are quoted at the "pivot" scale wavenumber of fc = 
0.05 Mpc-^ 

We assume that the reionization of the neutral inter- 
galactic medium occurred instantaneously, with an optical 
depth given by r. Finally, O gives the ratio of the sound hori- 
zon scale at the epoch of decoupling to the angular diameter 
distance to the corresponding redshift. 
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There are further basic quantities whose values can be 
derived from those listed in the set of Eq. ([H}: 

Pdorivcd = (r^m, /l, Ode, 0"8, 2rc, io)- (5) 

The matter density parameter is given by Q.m = Hdm + f^b- 
The value of the Hubble constant is derived from h — 
\/ (wdm + u!h)/i^in- The energy-density of the dark energy 
is set by Que = 1 — f2m — Qk- The results for As can be 
translated into a constraint on erg, the rms linear pertur- 
bation theory variance in spheres of radius 8 /i~^Mpc, us- 
ing the matter fluctuation transfer function. The redshift 
of reionization, Zra, can be computed from the values of r, 
the Hubble co nstant and the matter and baryon densities 
iTegmark et al. ..1994) . The age of the universe is to- 

The ACDM cosmological model is the simplest model 
which can account for the wide variety of cosmological ob- 
servations available today. This model is characterized by 
six parameters: 

Pvaricd = (t^dm,(^b,ns,T, ^s,e), (6) 

assuming f2k = and woe ~ —1- This parameter space is 
well constrained by the temperature and polarization power 
spectrum measur ements from flve years of integration of the 
WMAP satellite (|Komatsu et al. 1120091 '). In Section [STTI we 
analyse the impact of the measurement of the shape of ^(s) 

on the constraints in this parameter space. 

Using the latest WMAP data, iKomatsu et al. I (|2009l ') 
placed strong constraints on the possible deviations from 
the ACDM model, namely non-Gaussianity, the presence of 
isocurvature modes, deviations from a pure power law scalar 
primordial power spectrum, the presence of tensor modes, 
a non-negligible energy component in the form of massive 
neutrinos and parity- violation interactions. However, there 
are two parameters that signal important deviations of the 
standard ACDM model that can not be tightly constrained 
from CMB data alone: the curvature of the Universe Qk and 
the dark energy equation of state Wde- In order to assess 
the improvement on the constraints once the information 
on the shape of ^(s) is included in the analysis, we explore 
four parameter spaces which contain extensions of the simple 
ACDM set. First we extend the parameter set of Eq. (|6]) by 
adding a constant dark energy equation of state 

fixing fifc — 0. The results obtained in this case are shown 
in Section [5.21 

In Section FS.SI we include the parametrization of Eq. Q 
in our analysis and explore the extended parameter space 

Pviild"' = {iOdui,UJb,ns,T,As,Q,Wo,Wa), (8) 

where we also implement the hypothesis of a flat universe 
with f2fc = 0. We also analyse the effect of dropping this 
assumption. First, in Section 15.41 we include Qk as a free 
parameter, with 

Pvii^d = (i^dm,(^b,ns,r,yls,e,r2fc), (9) 

assuming that the dark energy is given by a cosmological 
constant (or vacuum energy) with woe ~ —1- Finally, in 
Section [5.51 we also allow this parameter to vary freely and 
we explore the full parameter space of Eq. ([3]). 



3.2 A physically motivated model for the 
correlation function 

In the linear perturbation theory regime (valid when the 
fluctuation amplitude is small, for instance at high redshift 
or on very large scales), the shape of the matter correlation 
function is well understood and can be readily obtained us- 
ing li near Boltzmann solvers, such as CAMB (|Lewis et al. I 
I2OO0I I. The shape of the correlation function, in the con- 
text of a standard adiabatic CDM model, is sensitive to 
the values of the matter density, Qmh/^, the baryon density 
O.bh'^, the spectral tilt ris and the density parameter of mas- 
sive neutrinos. The evolution with redshift of the correlation 
function is well understood in the linear regime. In this case, 
each Fourier mode of the density field evolves independently 
of the others and the shape of ^ is unaltered, although the 
overall amplitude changes with time. 

Unfortunately this simple behaviour is modified by a 
number of nonlinear phenomena, which affect different scales 
at different epochs. These include nonlinear effects gener- 
ated by the latter stages of gravitational instability, redshift- 
space distortions caused by gravitationally induced motions 
and a possible non-trivial scale dependent bias relation be- 
tween the distribution of galaxies and the underlying dark 
matter field (see, for example, the step-by -step illustration of 
these effects given bv lAngulo et al. II2OO8I ). Nonlinear growth 
results in cross-talk between different Fourier modes and in- 
troduces scale dependent patterns in the clustering, even on 
large scales. This is particularly noticeable in the appear- 
ance of the BAO bump, which is sensitive to the match be- 
tween the amplitude and phases of the fluctuations around 
the peak scale. Due to the distortion of the Fourier modes 
from their original values, nonlinear growth causes the BAO 
bump to become smeared out and also to lose contrast. Such 
effects cannot be ignored when modelling low-redshift data. 

On large scales, (r > 30 Mpc), the correlation func- 
tion falls sharply with increasing comoving pair separation, 
scaling roughly as a power law ~ r with 7 ~ 2.5. At even 
larger separations, this behaviour is altered by the emer- 
gence of a bump known as the BAO peak (e.g. see Fig. [3l). 
This f eature, first measured for LRGs by lEisenstein et al.l 
(|2005l ). is centered at about 110 Mpc and has a width of 
~ 20 Mpc. The amplitude of the peak in the LRG corre- 
lation function corresponds to a ~ 1% excess in the number 
of LRG pairs above the number expected in a random dis- 
tribution (e.g. Ca bre fc Gaztafi aga 2009a). 

To model the sh ape of the correlation f unctio n on 
large scales, we follo w ICrocce fc Scoccimarro I (120081 ) and 
I Sanchez et al. I l|2008l ) and adopt the following parametriza- 
tion: 

CNL(r) = &^{CL(r)®e-('=*'-)' +AmcCl(OC^,''M}, (10) 

where b, fc* and Amc are nuisance parameters, and the sym- 
bol (8) denotes a convolution. Here is the derivative of the 
linear correlation function and ^^\r) is defined by the in- 
tegral: 

C[''(r) = V"'eL(r) =47r J Pi^{k) ji{kr)kdk, (11) 

with ji{y) denoting the spherical Bessel function of or- 
der one. The model in Eq. [TD] is primarily motivated by 
RPT, where the matter power spectrum is written as P = 
G^Pl + PmCi with G a nonlinear growth factor and Pmc 
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Figure 3. The correlation function of dark matter halos in red- 
shift space measured in an ensemble of 50 large-volume N-body 
simulations (with total volume of ~ 105 Gpc^). The error 
bars correspond to the error on the mean of the ensemble and are 
obtained from the scatter among the 50 realizations. The best fit- 
ting parametric model used in this work, Eq. IIIOI I. is shown by the 
solid blue line (note the fit takes into account the covariance be- 
tween the bins). The dashed line corresponds to setting Amc = 
and highlights the importance of this term in matching the shape 
of the correlation function at r < 80 Mpc. The scaling uses 
r-BAO = 102/i-iMpc. 

being the power generated by mode-coupling. To a very 
good approximation G is of Gaussian form, while at large 
scales the leading order contribution of Pmc in real space 
is ~ Ci^'Ci (see Appendix |X] for a detailed discussion). 
The Gaussian degradation of the BAQ information was als o 
shown to be a good description bv lEisenstein et al. I l|2006h . 

We now demonstrate how accurately the model of 
Eq. pO|l can reproduce the spatial clustering of halo sam- 
ples with comoving number densities similar to that of the 
LRGs in DR6. To this end we utilize an ensemble of 50 re- 
alizations of coUisionless dark matter N-body simulations. 
Each simulation contains 640^ particles in a comoving vol- 
ume V = = (1280 /i"^ Mpc)^. The cosmological param- 
eters were set to Qm = 0.27, = 0.7, Qb = 0.046 and 
h = 0.72. The initial power spectrum had spectral index 
Us — 1 and was normalized to give erg = 0.9 when lin- 
early extrapolated to z = 0. Halos were identified using 
the friends-of-fri ends algorithm with linking-length param- 
eter / = 0.2 (see lSmith et al. Il2008l :lc rocce fc Scoccimarro I 
l2008l . for more details). 

Fig. |3] shows the 2-point correlation function measured 
in redshift-space ai z — for two non-overlapping samples 
of dark matter haloes of masses 7 X 10" < M[h-^MQ] < 
15x 10" and 15x lO" < M[h~^ Mq] (with number densities 
n[/i^Mpc"^] = 1.88 X lO"'^ and 3.46 x lO"'^ respectively). 
Our redshift space measurements were done considering the 
contribution of peculiar velocities along one dimension in the 
simulation (i.e. recreating a plane-parallel configuration). 



ISzapudi I (|2004l ) discussed in depth the issue of wide 
angle redshift space distortions. He found that opening 
angles 9 < 15° — 20° would ensure the validity of the 
plane-parallel approximation. For the redshift range of the 
LRG sample considered in our analysis (0.15 < z < 0.47) 
and the largest pair-distance separation allowed (smax = 
150 /I'^Mpc), the maximum opening angle is given by 

Sin(6lniax) = Smax/-DA(2 = 0.15), which implies 6'max « 20°. 

For the mean redshift of the sample the angle subtended by 
Smax is S « 10° , well withiu the validity of the plane-parallel 
approximation. 

The error bars plotted show the error on the mean ob- 
tained from the ensemble of 50 simulations. Each simulation 
has a larger volume than the LRG sample we use from DR6. 
Therefore these errors are « \/2 x \/50 — 10 times smaller 
than they would be for the LRG sample. This means that 
the deviations between models (solid lines) and simulation 
results (points) are not important for our purposes. 

The parametrization given in Eq. (|10|l corresponds to 
the solid blue line in Fig. [3] The corresponding best-fit 
values for the nuisance parameters (fc*[/i Mpc~^], &, j4mc) 
were (0.142,1.94,3.23) and (0.146,2.75,4.53) respectively. 
Clearly, the model can accurately describe the clustering 
of h alos in redshift s pace. 

I Sanchez et al. I (|2008l ) also compared Eq. (|10|l against 
measurements of the non-linear correlation function from 
a similar large ensemble of N-body simulations at various 
redshifts, and confirmed that this form gives an essentially 
unbiased measurement of the dark energy equation of state 
using both, real and redshift space information. 

The dashed line in Fig. [3] corresponds to Eq. (|10p with 
^MC = 0. From the plot we can see that the inclusion of this 
term is the key to recovering the correct clustering shape 
at separations smaller than the BAO bump. In addition it 
contributes slightly to the shape of the bump and alleviates 
a systematic effect related to the position of the BAO peak 
(|Grocce fc Scoccimarro |[2008l 'l (but this is sub-dominant for 
the survey volume being considered) . In subsequent sections 
we check that this nuisance parameter is not degenerate with 
any of the cosmological parameters. 

On theoretical grounds, one expects the smoothing 
length k^^ to depend on cosmology (i.e. aside from galaxy 
type or redshift) with, for example, a_10% increa se in Om ex- 
pected to increase fc^ by about 4% (jCroccc fc Scoccimarro I 
l2006bl : iMatsubara Il2008l ) . In view of this, we decided to con- 
sider fc+ as a nuisance parameter at the expense of a possible 
increase in error bars. Note that thi s is at variance w i th the 
approach of Eiscnstcin e t al. I l|2005h . pRgmark et al.1 (|2006l ) 
and Percival et al. (2007g) who kept this length fixed. 

In summary, for each cosmological model we compute 
the linear correlation function ^h{r) using CAMB to gener- 
ate the corresponding transfer function, and simulate non- 
linear effects through Eq. I|10p after computing ^L(r) and 
^i^\r) from Eq. (fTTIl. 



3.3 A model for the radial acoustic scale 

Here we describe the simple model that we use to compute 
the dimensionless redshift radial acoustic scale A^bao for a 
given choice of the cosmological parameters of Eq. . 
The value of Azbao can be computed as. 
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Table 1. The parameter space probed in our analysis. We assume 
a flat prior in eacii case. The parameter spaces that we consider 
are set out in Section |3!T1 



Parameter Allowed range 



e 



wde 



0.01 - 0.99 
0.005 - 0.1 
0.5 - 10 
- 0.8 
0.5 - 1.5 
2.7 - 4.0 
-2. - 
-0.3 - 0.3 



Azbao(-z) = 



(12) 



where H{z) is the Hubble constant at the mean redshift of 
the measurements (zm = 0.24 and 0.43), and rs(2d) is the 
comoving sound horizon at the drag epoch, which is given 
by 



da 



(13) 



where R = 'i^^/AVL^ and Vi^ = 2.469 x 10"^ /i"^ for a CMB 
temperature Tcmb = 2.725 K. The value of can be com- 
puted with high accuracy f rom the values of uih an d ojdm 
using the fitting formulae of lEisenstein fc Hul 1 19981 ). 

The radial BAO scale measurements can be used to 
place constraints on cosmological parameters independently 
of other datasets. Being a purely geometrical test, as in the 
case of the SN data, the radial BAO measurements are not 
sensitive to all the cosmological parameters of Eq. (|3]) since 
they contain no information about the primordial spectrum 
of density fluctuations, that is As and rig, and of the optical 
depth to the last scattering surface r. 



3.4 Practical issues when constraining parameters 

We use a Bayesian approach and explore the different pa- 
rameter spaces defined in Section 13.11 using the Markov 
Chain Monte Carlo (MCMC) technique. Our results were 
generated with the publicly available COSMOMC code of 
Lewis & Bridle (2002). CoSMOMC uses the CAMB package 
to compute power spectra for the CMB and matter fluctua- 
tions (Lewis, Challinor & Lasenby 2000). We use a general- 
ized version of CAMB which supports a time-dependent dark 
energy equation of state (^Fang ct al. 200^). For each pa- 
rameter set considered, we ran twelve separate chains which 
were stopped when the Gelman and Rubin (1992) criteria 
reached R < 1.02. We implemented flat priors on our base 
parameters. Table 1 summarizes the ranges considered for 
different cosmological parameters in the cases where their 
values are allowed to vary. 

In order to establish the link between a given cosmo- 
logical model and the datasets described in Section [2] it is 
necessary to include a small set of extra parameters given 
by 
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Figure 4. The impact of a mismatch in cosmology on the form 
of the correlation function. The upper panel shows the difference 
between the correlation functions measured assuming difli'erent 
fiducial cosmologies (ACDM models with varying Cm) from that 
obtained assuming Om = 0.25. The lower panel shows the same 
comparison after applying the scale shift of Eq. I I15I I to take into 
account the change in the value of Qm- 



(6, fc*,^MC,^Sz), 



(14) 



to the parameter sets described in Section 13.11 The bias 
factor b, describes the difference in amplitude between the 
galaxy correlation function and that of the underlying dark 
matter. The values of A;* and Amc from Eq. (I10|l are also 
included as free parameters in our parameter space. Asz 
gives the amplitude of the contribution from the Sunyaev- 
Zeldovich effect to the CMB angular power spectrum on 
small scales (high I). When quoting constraints on the pa- 
rameters of Eq. ([3} and ([5]), the values of these extra pa- 
rameters are marginalized over. In the case of b, this is done 
using the analytic expression given in Appendix F of Lewis 
& Bridle (2002), but the remaining parameters are included 
in our Monte Carlo analysis (see Section [3. 4p . 

There is one important point that must be considered 
in order to make a comparison between the model of Eq. (|10|l 
and the observational data of the LRG redshift-space cor- 
relation function. When measuring ^(s), in order to map 
the observed galaxy redshifts and angular positions on the 
sky into distances, it is necessary to assume a fiducial cos- 
mological model. This choice has an impact on the results 
obtained. 

One possible way to deal with this is to re-compute the 
correlation function and its covariance matrix for the cos- 
mology corresponding to each point in the Markov chains, 
and then use this measurement when computing the likeli- 
hood of the given cosmological model. This approach is in- 
feasible since it would require an exceedingly large amount of 
computing time. Instead we follow an alternative approach 
by indirectly including the effect of the choice of the cosmol- 
ogy on the model correlation function. To do this we follow 
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lEisenstein et al.l |200^ and simply rescale the distances in 
the model correlation function by a factor 



Dl 



^fiducial ' 



(15) 



where the effective distance Dy[z^) to the mean redshift of 
the survey z-m = 0.35, is computed for each model as 



-Dv(Zm) 



1/3 



(16) 



where Da{z) is the comoving angular diameter distance 
given by 

dz' 



Da{z) 



where 



fk{x) = 



/o H{z') 



smh{x) if(J7fc > 0), 
X if(!^fc=0), 
sin(a;) if(f2fe < 0). 



(17) 



(18) 



The exponents within the square bracket on the right- 
hand side of Eq. ((16} assume that the survey covers a wide 
solid angle, rather than a pencil-beam. 

We can test the effectiveness of the correction factor 
given by Eq. (I15|l by applying it to different estimates of the 
correlation functio n computed with varying cho ices of the 
fiducial cosmology. ICabre fc Gaztafiaga I ()2009ah computed 
the redshift space correlation function of SDSS-DR6 LRGs 
assuming different flat fiducial cosmologies with values of 
Qrn ranging from 0.2 to 0.3. The upper panel of Fig.|4]shows 
the difference of these estimates from the one obtained for 
Qui = 0.25. It is clear that the choice of the cosmological 
model affects the shape of the correlation function. The ef- 
fect of this choice is particularly important on small scales 
(r < 70 Mpc), but also introduces a small distortion 
to the shape of the acoustic peak, shifting its position to- 
wards larger (smaller) scales for smaller (larger) values of 
Qm- If unaccounted for, this difference would bias the con- 
straints obtained on the cosmological parameters. The lower 
panel shows the same quantity once the correction factor of 
Eq. (|15p has been applied. Clearly, this simple correction 
is able to account for the choice of the fiducial cosmolog- 
ical model. We use this correction factor to translate the 
mo del correlation function to t he fiducial cosmology used 
by ICabre fc Gaztanaea I ()2009ah to estimate the LRG ^(s) 
(a flat ACDM mode with = 0.25). We then compute 
the likelihood of the model assuming the Gaussian form 
£ cx exp(-xV2). 



4 TESTING THE MODEL OF C(S) 

In this section we analyse the sensitivity of the parameter 
constraints to the details of the procedure we follow to com- 
pute the likelihood of a given model. For the purpose of this 
exercise we use the information contained in the correla- 
tion function of LRGs combined with the latest results from 
the WMAP satellite alone to constrain the parameter set 
of Eq. ([7|, and assess the impact on the results of varying 
choices in our analysis procedure. 

First we study the sensitivity of the resu lts to the range 
of scales in ^(s) included in the analysis. ISanchez et al. I 
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Figure 5. The dependence of the constraints on cosmological 
parameters on the minimum pair separation, rmin, included in 
the correlation function measurement. The points show the mean 
value of the likelihood for each parameter and the error bars show 
the 68% CL. 



(l2008t ) showed that the model of Eq. HOJ gives an accu- 
rate description of the redshift-space halo correlation func- 
tion measured from N-body simulations on scales in the 
range 60/i"^Mpc < r < 180/i"^Mpc. In Section [32] we 
showed that the second term in Eq. (|10p also helps to re- 
produce the results from N-body simulations down to scales 
of 40/i~^Mpc. Even though our measurement of ^(s) ex- 
tends to larger scales, the measurement beyond 150 Mpc 
is noisy and does not have any effe ct on the constraints. 
Besides, ICabre fc Gaztafiaga I (j2009ah showed that on these 
large scales, the measurement of the redshift space correla- 
tion function exhibits an excess in the a mplitude with re - 
spect to the previous estimate of Eisenst ein et al. I (|2005l ). 
For this reason we have chosen to set a maximum distance 
^max = 150 Mpc in all subsequent analyses. 

Fig. [5] shows the one-dimensional marginalized con- 
straints on a subset of the parameters space of Eq. ([7]) as a 
function of the minimum scale, r-min, included in the analy- 
sis. For r-min < 40 Mpc we see that the constraints start 
to vary with rmin. Coincidentally, in this regime the model 
of Eq. (I10|) starts to deviate from the N-body measurements 
of the correlation function discussed in Sec. 13.21 On such 
scales, extra terms beyond the one-loop contribution ~ ^' 
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Table 2. The marginalized 68% interval constraints on the cosmological parameters of the ACDM model obtained using different 
combinations of the datasets described in Section |2] as stated in the column headings. 





CMB 


CMB + 


CMB + rBAO 


CMB + SN 


+ rBAO 


owirs + ^ i^sj 
+ rBAO + SN 


rBAO 


looe 


1.0416+o»»23 

— 0.0023 


1.0415+Hn22 

— 0.U022 


1.0413+nn22 

— O.UOzz 


1.0413+o°n?^ 


1.0414+H!!^o 

— O.UOzz 


1.0412+O-On?? 

— U.0022 


0.987+0-060 

— 0.069 




0.1091+°°»^3 

^ 0.0053 


0.1076+°°°38 

v^. j-v< 1 '-'—O.OOSS 


0.1108+°0036 

^ ^^—0.0036 


n ii2c;+o oo3s 

^ '-^— 0.0038 


0.1095+°0032 

"^'^^— 0.0032 


0.1110+00029 

^ ^— 0.0029 


0.091+0-040 

^^'^— 0.038 


100 uJh 


■ —0.051 


2.280+Hm 

— 0.051 


2 277^^-^^^ 
— 0.050 


2.270+0°^l 

— 0.051 


2.276+°°4? 
— o.04y 


2.267+0-050 

— 0.051 


3.0+?-0 

— 1.7 


T 


0.088+°»i^ 

"^^^—0.016 


0.090+0°;^ 

— 0.017 


0.088+°°|8 

^^^^—0.016 


0.087+?,-i!l^ 

— 0.016 


0.088+°°l« 

—0.016 


0.086+0 016 

^^^"—0.016 






0.9651H13 


0.965lHi2 


0-962lHi2 


0-9601H12 


0-963lHn 


o-seolHll 


i-06lS;^t 


In(lOlOAs) 


J.UDD_Q 


3.059lg°f, 


3 067+" °3'^ 


3 070+° "^'* 

O.UlU_Q g3g 


3.062l°:°36 


3 065+0 036 
O.UDO_o Q3g 






n 740+0.026 


n 7Cip;+0.018 


n 740+0 017 


U./OU_Q 020 


n 747+0.014 
'^'-0.014 


(, -oq+0.013 
u./jy_o 0^3 


n 709+0. 040 




„ r,c:i +0.026 

U.^O±_Q Q2g 


n 244+0 018 
^•^^^-0.018 


0.260lHi^ 


27n+° °2o 




0.26110:013 


n 91Q+0.041 


erg 


0.795lH^^ 


0.787lH^2l 


0.802lH^,l 


o.siolHS 




0.80210:0^1 




to/Gyr 


13.64t«- 


13.6410:1° 


13.668l°:0f, 


13.69l°:l« 


13.66010:0^^ 


13.68210:095 


14 4+3-2 




10.511:^ 


10.511:3 


10.411:^ 


10.411:3 


10.411:3 


10.411:3 




h 


0.726l°;[Jg 


r, 701+0.018 

U. (01_Q QJg 


0.7181^^6 


n 7in+0-0l7 


f^ 790 + 0.013 


0.71610:012 


n 70+0.15 
''^-0.14 



to ^MC are expected to become important, which implies 
the breakdown of our ansatz. 

For rmin > 40/i~^Mpc, Fig. [S] shows that the con- 
straints on the values of the cosmological parameters are 
very stable. Furthermore, the mean values obtained for these 
parameters using CMB data plus the LRG correlation func- 
tion are in complete agreement with the ones obtained from 
CMB information only. The allowed regions of some parame- 
ters, like tJb or rig which are tightly constrained by the CMB 
data alone show almost no change on varying rmin. Other 
parameters, such as Sim or iude show a substantial increase 
in their allowed regions as the data from small scales is grad- 
ually excluded from the analysis. For rmin ~ 42.5 Mpc 
we obtain a constraint on the dark energy equation of state 
of WDE = -0.996l!l:!1^5, while for rmin = 82.5 h'^ Mpc, that 
is including only scales close to the acoustic peak in ^(s), 
we get Woe = — CQQlg lg. This approximately factor of two 
change in the error bar highlights the importance of the 
inclusion of information from the shape of ^(s) on interme- 
diate scales. Based on this comparison and the results of 
Section [3.21 we have chosen to set rmin ~ 42.5 Mpc in 
the subsequent analysis. 

We also studied the impact of the correction to the mea- 
sured correlation function to take into account changing the 
fiducial cosmology assumed when estimating ^(s). For this 
test we obtained constraints in the same parameter space 
without applying the correction of Eq. (|15p . The results ob- 
tained in this way are entirely consistent with those obtained 
when this correction is applied. The mean values of all the 
parameters remain practically identical with the exception 
of Qm which shows a slight shift towards higher values of 
approximately 0.2a. The allowed regions for these parame- 
ters also show almost no variation with a slight increase in 
the confidence limits of Qm and wbe of approximately 0.3(7. 
A similar shift but towards smaller values of ilm is obtained 
on setting Amc = in the second term of the right hand 
side of Eq. ([101). Th is is in agreement with the results of 
ISanchez et al. I l |2008l ). who found that the two most impor- 



tant parameters of this parametrization required to obtain 
a good description of the shape of the two-point correlation 
function are fc* and b (which in this case is marginalized 
over) . 

In Section [S] we also present constraints obtained from 
the combination of the two large scale structure datasets we 
use in our analysis, the LRG ^(s) and the position of the 
radial acoustic peak, without including any CMB informa- 
tion. The fact that the correction for the fiducial cosmology 
does not significantly change the obtained results when the 
LRG 5(s) is used in combination with the CMB, is not due 
to the later keeping the constraints sufficiently close to the 
fiducial cosmology. As we will see later (see Section [5. 2|l the 
constraints in this parameter space from CMB data alone 
present a strong degeneracy between fim and tode, allowing 
for models that differ substantially from our fiducial cosmol- 
ogy of flm ~ 0.25 and wue = — 1- Nonetheless, these models 
are ruled out by the ^(s) data itself, (which on large scales 
is insensitive to this correction), since the position of the 
acoustic peak shows a strong variation in these models. The 
only region of the parameter space allowed by the data is 
that where the effect of the correction of Eq. p5p is small 
(although not completely negligible) . This is the reason why 
the obtained constrains are not extremely sensitive to this 
correction. This also implies that this correction should not 
be too important even in the case in which the LRG ^(s) 
is combined with the radial acoustic peak whitout including 
any CMB data. We have tested this explicitly and found 
that the resulting mean values of the cosmological param- 
eters obtained for this dataset combination when this cor- 
rection is not applied show no variation with respect to the 
ones obtained when it is used. The only effect of ignoring 
this correction is a small artificial increase of the allowed 
regions of Sim and wue of the same order as in the previous 
case. Setting Amc = for this combination of datasets also 
produces a small shift towards smaller values of Sim, which 
in this case is much smaller than the allowed region for this 
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Figure 6. The marginalized, one-dimensional posterior likelihood in the ACDM parameter space (Eq.O obtained from CMB information 
only (dashed lines), CMB plus the shape of ^(s) (solid line) and the full constraints including also rBAO and SN data (dot-dashed lines). 



parameter. These tests show that our constraints are robust 
with respect to the details of our analysis technique. 



5 CONSTRAINTS ON COSMOLOGICAL 
PARAMETERS 

In this section, we carry out a systematic study of the con- 
straints placed on the values of the cosmological parameters 
for the different parameter spaces defined in Section [3.41 In 
Section 15.11 we present the results for the simple ACDM 
cosmological model with six free parameters. In Section [5.21 
we consider an extension to this parameter set, allowing the 
dark energy equation of state iude to float (but without any 
redshift dependence) . Section 15.31 gives the constraints on 
models where the time variation of wde is parametrized ac- 
cording to Eq. (|4]). In Section [5]4] we discuss our constraints 
on non-fiat models, analysing the parameter space of Eq. Q. 
Finally, Section [5.51 shows the constraints on the full param- 
eter space of Eq. ([3)l , allowing both for non-flat models and 
more general dark energy models. Tables I2I6I compare the 
constraints obtained in these parameter spaces using differ- 
ent combinations of the datasets described in Section |2] 



5.1 The basic ACDM model 

Due to the successful reproduction of a wide variety of obser- 
vations, the ACDM model has emerged over the past decade 
as the new standard cosmological model. The recent results 
from five years of observations by the WMAP satellite have 
helped to reinforce this conclusion. The latest WMAP data 
give a much better estimation of the third acoustic peak in 
the CMB temperature power spectrum, as well as the low- 
l polarization signal. Thanks to these improvements, the 
WMAP data alone has been able to provide much tighter 
constraints on this basic cosmological model than was pos- 
sible when using earlier releases. 

Table [2] summarizes the constraints on the parameters 
of this simple model for different combinations of datasets. 
Fig. [6] shows the marginalized likelihoods for this param- 
eter set obtained using CMB data alone (dashed lines), 
CMB plus the LRG ^(s) (solid lines) and the full combina- 
tion of the datasets described in Section [2] i.e. CMB-I- LRG 
^(s)-l-rBAO-l-SN (dot-dashed line). Several parameters, such 
as ti;b, tJc and r, are tightly constrained by the CMB data 
alone, and show almost no variation when other datasets are 
included in the analysis. On the other hand, the constraints 
on the parameters of the energy budget do show a marked 
improvement on adding further datasets. 

The CMB data is particularly sensitive to the value of 
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Table 3. The marginalized 68% interval constraints on cosmological parameters allowing for variations in the (redshift independent) 
dark energy equation of state (i.e. the parameter set defined by Eq.^Tj), obtained using different combinations of the datasets described 
in Section [21 as stated in the column headings. 





CMB 


CMB + F(f) 


CMB + rBAO 


CMB + SN 


CMB + 
+ rBAO 


CMB + 5(s) 
+ rBAO + SN 


rBAO 


Woe 


n 70+0.30 

'^■'■^-0.30 


n QSS+O-088 


-0.92t0-;j 


-0 950+0-054 
U.yOU_o 055 


-0 999+0-090 
u.ayj_o 091 


-0 969+0-052 
u.yDy_o 052 


-1 05+0-16 
l.'-"^-0.15 


looe 


-10410+0.0023 

iU41J_Q QQ24 


1 0415+0 0023 
i.U4iO_o 0023 


l-0414t»:«i 


1 0413+0-0023 
i.U41J_o 0022 


1 0413+0-0023 
i.U41J_o 0023 


1.041410.0022 


994+0.071 
U.Ma4_o 075 






n infi8+o oo49 

U.lUfaS_o.oo48 


1095+0-0045 
U.iuyo_o 0045 


10033+0.0049 
u.iuyoo_o 0050 


1092+0-0042 
u.luyz_o 0042 


0.1088lg:°°l? 


099+0.049 
u.uya_o 045 


100 OJ^, 




9 9on+0 052 
^.^»U_o 052 


9 977+0.051 
'-0.051 


9 979 + 0.052 
Z.Z(Z_o 052 


9 97C-+0. 00051 
Z.Z(0_o 049 


9 97C+0.051 
z.Z(0_o 050 




T 


0.089to.oir 


090+0 017 
u.uyu_o 017 


090+0-017 
u.uyu_o 017 


o.osolHl^ 


0.088t0-0i7 


0.088lHi«6 


- 


ns 


959+0 014 
u.9oy_g Qj^4 


965+0-012 
u.yDO_o 013 


963+0-012 
u.»uo_o 012 


0.9611^^3 


0.963««li 


963+0-012 
u.yoo_o 012 


1 03+0.37 
l.^'-'-0.32 


In(lOiOAs) 


067+0-038 


nt^7+0.037 
o.uol _o.o37 


3 066+0-037 
J.UDD_o 037 


3 064+0-038 
O.UO'i_o 037 


3 062+0-035 
o.uoz_o 036 


3 060+0-036 
d.UDU_o 036 






n 00+0.12 

0.631,) 12 


n 7^4+0.020 

U. (04_o 021 


PI -90+0.034 
U. '^o_o 034 


n 700+0.018 
U. 'oo_0.019 


746+0-017 
U./4D_o 017 


700+0.013 

u. My_o 013 


n 770+0.045 
0-/ 'S_o.045 




S6+012 


n 245+0-021 
u.zio_o 020 


n 9-7+0.034 
"J-^' ' -0.034 


0.267lHf8 


254+0-017 
U.Z04_o 017 


261+0 013 
U.ZDi_o 013 


n 999+0.045 
U.ZZZ_o 045 




'J- '^^-0.085 


0.778l°:°li 


n 774+0.060 

'J- ' '^-0.060 


0.7801^^7 


703+0.043 
U./yj_o 045 


0.78ltHl! 




to/Gyr 




13.65«:11 


13 74+0-14 
l'^- '^-0.14 


13.7lt«:i° 


13.67lg:ig 


13.689to«^^ 


14 3+^-^ 








in fi + 1-4 
10.6^1.3 






io.4li:^3 




h 






695+0-046 

U.DMO_o 046 


0.704lo.or7 


r, 799+0.024 
'^^-0.025 




74+0.16 
'J- '*-0.15 



uJm = flmh^- This leads to a degeneracy between Qui and h column of Table [2] In this case we get f2n 



: 0.218l° «^J and 



which can be broken b y combining the CMB measurements 
with other datasets (|Percival et al. I l2002l: ISoergel et alTI 



h = 0.73: 



Although these constraints are weaker than 



I2OO3I : ISanchez et al. I I2OO6I : ISpergel et al. II2007I ). However, 
the improved estimation of the third acoustic peak in the 
WMAP5 temperature power spectrum helps to alleviate the 
degeneracy with respect to earlier data releases. The ratio 
of the amplitudes of the first and third acoustic peaks in 
the temperature power spectrum is sensitive to the ratio 
Qui/^i- This results in an improvement in the constraints 
on Qui from CMB data alone, thereby reducing the degen- 
eracy between Qm and h and improving the constraints 
on these parameters. From the CMB data alone we get 
fim = 0.251 ± 0.026 and h = 0.726t°o°oll- Including the 
shape of C(s), these constraints change to Qm = 0.244±0.018 
and h = 0.731 ± 0.018, in complete agreement with the re- 
sults from the CMB and with previous determinations based 
on the combina,tion o f CMB and large scale structure data 
/Sanche z et al. I I2O O6': 'Spergel e t 811120071 : iDunklev et al. I 
[2009. : .Komatsu et al. 2009). The data from the radial BAG 
and SN prefer slightly higher values of Qm than the CMB 
data and the LRG ^{s), but are consistent within la. Com- 
bining the information from all these datasets we get our 
tightest constraints, with Qm — 0.261 ± 0.013 and h = 
0.716 ±0.012. 

With the exception of the optical depth r and the am- 
plitude of density fluctuations As, it is also possible to ob- 
tain constraints on the same set of cosmological parame- 
ters from the combination of the LRG ^(s) with the rBAG 
data without including any CMB information. For this we 
have included the bias parameter b explicitly in the anal- 
ysis, but marginalized the results over a wide prior with 
0.5 < b < 20. This prior has minimal impact on the obtained 
constraints but it allowed us to effectively test if a prior in b 
can have implications on the obtained constraints. The re- 
sults with this combination of datasets are shown in the last 



those obtained using CMB data, their importance lies in the 
fact that they are determined purely on the basis of large 
scale structure information. 

Gn co mbining the WMAP da. ta with the BAG measure- 
ment s fromlPercival et al. I |2007d) and the same SN UNIGN 
data. lKomatsu et al. I(|2009l ) found Qm = 0.279±0.015. This 
value is only marginally consistent with our results for CMB 
plus 5(s). This might indicate systematic problems intro- 
duced by the approximate treatment of the BAG measure- 
ments in previous analyses. We shall return to this point in 
Section [T] 

From the analysis of a compilation of CMB mea- 
surcmcnt s with the final power spectrum of the 2dFGRS, 
ISanchez et al. I l|2006il found evidence for a departure from 
the scale invariant primordial power spectrum of scalar fluc- 
tuations, with the value Us = 1 formally excluded at the 95% 
level. This deviation was subsequently confirmed with higher 
significance wi th the avai lability of the three-year WMAP 
data (Spergel et al. I [20071 ). By the full combination of the 
datasets of Section[2]we find a constraint on the scalar spec- 
tral index of = 0.963to.m}. 

with the Harrison-Zel'dovich 
spectrum 3.6cr away from the mean of the distribution. 

The conclusion from this section is that the ACDM 
model gives a consistent and adequate description of all the 
datasets that we have included in our analysis. The precision 
and consistency of the constraints on the basic parameters 
on this model constitute a reassuring validation of the cos- 
mological paradigm. In the following sections we will concen- 
trate on two possible deviations from this model that can be 
better constrained by the shape of £,{s), namely alternative 
dark energy models and non-fiat cosmologies. 
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Figure 7. Panel a): the marginalized posterior likelihood in the 
f2in — «)DE plane for the ACDM parameter set expanded by the 
addition of iode (Eq. 10). The short-dashed lines show the 68 
and 95 per cent contours obtained using CMB information alone, 
solid contours show CMB plus LRG ^(s) constraints. Panel b) 
Comparison of the marginalized posterior likelihood in the same 
parameter space obtained using CMB information plus LRG 
(solid lines), CMB plus the radial BAO signal (long-dashed lines) 
and CMB-I-SN (dot-dashed lines), as indicated by the key. The 
filled contours correspond to the 68% CL in each case. 



5.2 The dark energy equation of state 

When treated as standard candles, the apparent dimming 
of distant Type la supernovae surprisinRl y pointed towards 
an accelerating expansion of the Universe (jRiess et al. II 199^ : 



xp, 

IPerlmutter et al. ll 999':'Rie ss et al. |2004| ). This was the first 
piece of observational evidence in favour of the presence of 



Figure 8. The marginalized posterior likelihood for the dark en- 
ergy equation of state parameter, judEi in th^ case of the parame- 
ter set defined by Eq. lO . The dashed line shows the likelihood in 
the case of CMB data alone, the solid line shows the likelihood for 
CMB combined with the LRG ^(s) and the dot-dashed line shows 
the likelihood using all of the data sets described in Section 2. 



a negative pressure component in the energy budget of the 
Universe. Independent support for this component, called 
dark energy, came from the combination o f CMB measure- 
ments and large scale structure data (Efstathiou et al. I2OO2I : 
iTegmark et al. II2OO4I ). Understanding the nature of dark en- 
ergy has become one of the most important problems in 
physics today since it has strong implications for our under- 
standing of the fundamental physical laws of the Universe. 
The simplest possibility is that the dark energy corresponds 
to the vacuum energy, in which case it behaves analogously 
to Einstein's cosmological constant with wde = —1, but 
several alternative models have been proposed. One way to 
narrow down the wide range of possible models is to obtain 
constraints on the dark energy equation of state parame- 
ter uiDE- In this section we extend the parameter space of 
the ACDM models to allow for variations in the (redshift- 
independent) value of uide. In Section [5.31 we drop this hy- 
pothesis to explore the possible redshift dependence of this 
parameter. Table |3] summarizes our constraints on this pa- 
rameter set from different combinations of the datasets de- 
scribed in Section [2] 

Fig. [TK shows the two-dimensional marginalized con- 
straints in the ilm— wde plane from CMB data alone (dashed 
lines) and CMB plus the LRG ^(s). There is a strong de- 
generacy between these parameters when only CMB data is 
included in the analysis which leads to poor one-dimensional 
marginalized constraints of wde = — 0.731q3q and Vlu-i ~ 
O.SGIq 5^2- Another view of this is given by Fig. [S] which 
shows the one-dimensional marginalized constraint on wde- 
The CMB only case is again shown by the dashed line. 

The origin of this degeneracy is well understood. The 
position of the acoustic peaks in the CMB power spectrum 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological parameters from the SDSS-DR6 LRGs ^(s) 15 



depends on the size of the sound horizon at the decoupling 
epoch, rs{Zf), which is given by Eq. (|13|l . The mapping of 
the physical scales of the acoustic peaks to angular scales 
on the sky depends on the comoving angular diameter dis- 
tance, Da (2*), given by Eq. (fT7)) . Therefore the peak pat- 
tern in the CMB pr ovides tight const r aints on the "acous- 
tic scale" given by feond et al. Ill99f : Efstat hiou fc Bond I 
1 19991 : iPage et al. H2003I : iKomatsu et al. „200a ) 



£a 



7vDa{z* 



(19) 



While wbe is relevant for the calculation of Da{z,), it has 
minimum impact on rs{z,), since the dark energy is dynam- 
ically negligible at decoupling. For this reason, for fixed val- 
ues of LJh and Wdm, and given a value of (or h), it is 
always possible to find a value of woe such that the value 
of £a remains constant. This gives rise to the degeneracy 
between these parameters seen in Fig. [7^. 

As shown in Fig. [7^, the inclusion of the LRG correla- 
tion function breaks the degeneracy between Sim and woe 
present in the CMB data. This is done in two ways; first the 
shape of ^ (r) on intermediate scales tightens the constraints 
on ilm, which helps to break the degeneracy between this 
parameter and woe- Second, through the position of the 
acoustic peak it provides an independent estimation of the 
ratio rs{zd)/Dv{zni), where z^ is the redshift of the drag 
epoch and 2m = 0.35 is the mean redshift of the survey. 

When this information is combined with the constraint 
on rs(zd) provided by the CMB data, this provides an extra 
distance measurement, _Dv(zm = 0.35), which breaks the 
degeneracy in the CMB data. This can be seen more clearly 
in Fig.|5J which shows the two dimensional constraints in the 
plane woE — Dv{zin = 0.35). Varying woe to keep a constant 
value of £a produces varying values of Dvizm = 0.35). The 
extra information from the shape of ^(s) fixes the value of 
Dv, tightening the constraints on the dark energy equation 
of state. In this case we get Qm = 0.245lo 020 ^'^'^ "Wde = 
— 0.9881o'q88' complete agreement with the cosmological 
constant. Again, Fig. |8] shows the dramatic reduction in the 
width of the likelihood distribution for wbe on combining 
the CMB data with the measurement of the LRG correlation 
function. 

Using WMAP data combined with the BAO m e asure - 
ment from iPercival et al. I l|2007d ). iKomatsu et al. I (|2009l ') 
found wde ~ — l-15to'22 (68% CL). If we exclude the small 
scale CMB experiments and consider only WMAP measure- 
ments plus the LRG £,{s), we get wde ~ — 0.996j;[5^Qg5, 
which corresponds to a reduction of almost a factor two 
in the allowed region for the equation of state parameter. 
This highlights the importance of the information contained 
in the shape of the c orrelation function. The method of 
iPercival et aTl (|2007d ) sacrifices the long-wavelength shape 
of P{k), which is affected by scale-dependent effects, in or- 
der to obtain a purely geometrical test from the BAO oscil- 
lations. 

It is also possible to obtain constraints in this param- 
eter space from the combination of the shape of the LRG 
^(s) with the rBAO data without including any CMB in- 
formation. In this case we get we get wde = — 1.051° ;[g. 
This means that present day observations allow to obtain a 
competitive constraint on iude purely on the basis of large 
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Figure 9. The marginalized posterior likelihood in the D^j —wj^^ 
plane for the ACDM parameter set extended to allow Wy>e to vary 
(Eq. (O). The dashed lines show the 68 and 95 per cent contours 
obtained using CMB information alone. The solid contours corre- 
spond to the results obtained from the combination of CMB data 
plus the shape of the LRG ^(s). The filled contours correspond 
to the 68% CL in each case. 



scale structure information, independently of CMB or SN 
observations. 

The values of uide listed in Table[3]as being obtained by 
combining the CMB data with rBAO and SN information 
are completely consistent (see Section [7] and Fig. [jls). Once 
again this shows the consistency between these datasets. 
Our results includin g rBAO data show an exc ellent agree- 
ment with those of ICaztafiaga et al71 l|2008 j ). When in- 
cluding systematic effects in the SN data as advocated by 
Kowalski et al. (2008), the CMB plus SN result changes to 



W-DE 



-0.900 



+0.078 



, a 50% increase in the allowed region 
and a shift in the mean value of wdb of about one a. This 
means that when the SN systematic errors are included in 
the analysis, the precision of the constraint on wde obtained 
from CMB plus SN is comparable to the one derived from 
CMB plus the shape of the LRG C(s). This highlights the 
importance of using the full shape of ^(s) as a cosmological 
probe. This also shows the importance of a precise determi- 
nation of the effects of systematic errors in the SN data. 



5.3 The time evolution of iude 

In the previous section we analysed the possibility of extend- 
ing the ACDM model with alternative dark energy mod- 
els with a redshift-independent equation of state. From a 
theoretical perspective, if wde 7^ —1, there is no real rea- 
son why it should be constant. In this section, we analyse 
the constraints on the redshift dependence of this param- 
eter, using the popular linear parametrization of Eq. Q 
JChevallier fc Polarski ir200ll : iLinder et al~ll2003l ) . 

This parametrization may lead to models in which the 
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Figure 10. The marginalized posterior likelihood in the wq — Wa 
plane for the ACDM parameter set extended with a redshift de- 
pendent dark energy equation of state parametrized according to 
Eq. l(4]l. The dashed lines show the 68 and 95 per cent contours 
obtained using CMB information alone. The solid contours cor- 
respond to the results obtained from the combination of CMB 
data plus the shape of the LRG S,{s). The dot-dashed contours 
show the constraints on using all datasets. The filled contours 
correspond to the 68% CL in each case. The dotted straight line 
marks the degeneracy between too a^nd ^a- 



Figure 11. The marginalized posterior likelihood in uide ^ ^ 
function of z for the ACDM parameter set extended with a red- 
shift dependent dark energy equation of state parametrized as in 
Eq. Q . The dashed lines show the mean value (thick line) and the 
corresponding 68 per cent constraints (filled blue region between 
thin lines) obtained using CMB information alone. The solid lines 
(and filled orange region) show the corresponding results from the 
combination of CMB plus the shape of the LRG J(s). The dotted 
vertical line corresponds to the pivot redshift 2p = 0.39. 



dark energy density has a dynamical impact at the epoch 
of Big Bang Nucleosynthesis (BBN), unlike the case with 
a cosmological constant. Such a scenario would affect the 
present day abundances of light elements, thereby violating 
the constraints on the possible variation in the Hubble pa- 
rameter at the ep och of BBN as expressed through the ratio 
l|Steigman 1120071 '): 



-ff'(aBBN) 



0.942 ± 0.030, 



(20) 



where -ff(aBBN) is the Hubble parameter in a standard 
ACDM model, in which case dark energy is completely neg- 
ligible at the expansion factor corresponding to the epoch of 
Big Bang nucleosynthesis, obbn, and H' [a^w) is the Hub- 
ble parameter at the same epoch in an alternative model 
in which dark ener gy do es pla y a r ole at this early epoch . 
Here we follow iWright] (l2007t ) and iKomatsu et al. I (120091 ') 
and impose a Gaussian prior given by Eq. (|20|) on 5* which 
can be written as 



S- 



O „-3(l+™off (zbbn)) 



(21) 



where abbn = 10"^ and WcB is the effective dark energy 
equation of state, defined as 



1 



ln(a) 



WDE(a') d In a . 



ln(a) Jo 

With the parametrization of Eq. (|4]) this becomes 



(22) 



lOcff (a) = wo + Wa { I + 



ln(a) 



(23) 



Table |4] summarizes the constraints obtained in this param- 
eter space using different combinations of datasets. 

The dashed lines in Fig. [10] show the two dimensional 
marginalized constraints in the wo — Wa plane from CMB 
data alone. As discussed in the previous section, the CMB 
data follows a degeneracy of constant ^a- Allowing the dark 
energy to evolve with redshift, this degeneracy gains an extra 
degree of freedom leading to poor constraints on wo and Wa ■ 

These constraints can be transformed into constraints 
on the value of wo^{a) which will be given by 



{5wu¥.{a) ) = {{5wo + (1 - a)5w^) ). 



(24) 



The dashed lines in Fig. [TT] show the mean value of 
w-d-e{z) as a function of z (thick line), as well as the 68 
per cent confidence limits (thin lines) obtained from CMB 
data alone. These constraints can accommodate large varia- 
tions of u)de(z). The solid lines show the correspondent re- 
sults when the LRG correlation function is included in the 
analysis. Although considerable deviations from the simple 
cosmological constant case are still allowed, these are much 
strongly constrained by the data. 

The pivot scale fa ctor Op is defined a s the point where 
Eg. (I24l) is minimized (iHuterer fc Turner ll200I : .Hu fc Jain I 
l2004l : lAlbretch et al. II2006I ). That is 



1 + 



(Swi) ■ 



(25) 
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Table 4. The marginalized 68% interval constraints on cosmological parameters allowing for an evolving dark energy equation of 
state (i.e. the parameter set defined by Eq.l[8]l), obtained using different combinations of the datasets described in Section[2l as stated 
in the column headings. 





CMB 




CMB + rBAO 


CMB + SN 


CMB + 
+ rBAO 


CMB + ^{s) 
+ rBAO + SN 


ffs'ia. rBAO 


wo 


-o.est-ll 


-l.OOiO-Jo 


-n 85+0-49 
U.OO_o.49 


-1 07+0' 10 
l-'J' -0.15 


-1 10+0.10 


-1.0310:10 


-1 Ofi+o.io 

l.'J"-0.18 




-0.30t«;«« 




-0 2+10 


45+0-59 


n 01 +0.60 
'.'■•^1-0.59 


0.3010:11 


-o.oin 


looe 


1 0413+00023 

l.U'HO_Q QQ23 


1 0416+0 0023 


1 0415+0 0022 

l.U'HO_Q QQ22 


1 0413+0 0022 

l.U'HO_Q QQ22 


1 0414+0.0022 
l.U'ilt_o 0022 


1 0416+0 0022 

l.U41D_o 0022 


997+0.073 
u.yy / _o 078 




o.ioy8_p 0054 


1076+0 0051 


n inQS+0 0043 

0.1098_o P044 


1107+0 0053 


n in«7+o.oo47 

O.i0»l'_o 0047 


0.107710:0^4^ 


096+0 048 

U.U9D_o 045 


100 Uli, 


9 97n+0-051 


2 284+0-053 
z.zot_Q Q52 


9 97Q+0.051 
^•^'"-0.051 


9 971 +0.053 
^■^'^-0.053 


9 97Q+0.051 
z.zio_o 051 


2.28010:01? 


2 e+i-** 

^-°-1.5 


T 


ogo+0017 


090+0 017 


090+0 017 
u.uyu_Q Qj^y 


0.08910-01? 


0.08910:01? 


091+0 016 
u.uyi_o 017 


- 


Us 


961+0-013 


966+0 013 


0.96410.012 


0.96010:0^3 


o.oeslHl^ 


963+0.012 
u.»uo_o 012 


1.0410:32 


ln(10i°As) 


3.0651003B 


3 o^q+0 038 


3.068l0;03| 


3.06810 037 


3 Ofil +0.037 

o.uui_o 037 


3.06210:036 








n 7c:c+0.023 
u. 1 JJ_o.023 


U. (10_Q Qg2 


n 718+0.024 
U-/18_o.024 


749+0.018 
U. '49_o 018 


0.73810:0" 


0.78010:°!^ 




34+0 " 


245+0 023 


285+0 062 

U.ZOO_g Qg2 


n 9C9+0.024 
u.zoz_o 024 


251+0.018 
U.zoi_o 018 


262+0.013 
u.zoz_o 013 


220+0-047 
U.ZZU_o 045 


OS 


p. „oo+0.084 

U. /00_Q gg2 


0.786l0;0« 


n 770+0.071 

'J-O.OTO 


7790+0 041 


791+0.047 
u. /3l_o 048 


769+0.040 

U.(D9_o 040 




io/Gyr 


13.89t0;32 


13.65tO;« 


13.72to.i4 


I3.83l0;l® 


1 79+0. 13 
'^-0.13 


13.7610:1^ 


14 6+3-3 






io.5t;* 


I0.5t;;l 


10.611:1 


10.511:1 


10.611:3 




h 


0.65tO;io 


r, 7qr, + 0.028 

U. (OZ_g 


0.6921^^5 


0.68910 023 


0.726l0:0?7 


706+0-015 

U./UD_o 015 


70+0. 16 
'■^-0.14 



For the combination of CMB plus the shape of ^(s) the cor- 
responding pivot redshift is given by Zp = 0.39, -which is 
sho-wn by the dotted line in Fig. 1111 At this redshift we get 
our tightest constraint on the dark energy equation of state, 
with wde(2p = 0.39) = —1.01 ± 0.10, entirely consistent 
with a cosmological constant. In this case also, adding the 
information from the shape of ^(s) allows to obtain a tight 
constraint on fim = 0.245lo:o22- 

The solid lines in Fig. [TO] show the two dimensional 
marginalized constraints in the wo — Wa plane for CMB plus 
the shape of ^(s). The constraint on wde(2p = 0.39) corre- 
sponds to a degeneracy between these parameters approxi- 
mately given by wq + 0.28wa = — 1, which is shown by the 
dotted line in Fig. \W\ 

Adding information from rBAO or SN gives completely 
consistent results, which are shown by the dot-dashed lines 
in Fig. 1101 The combination of all the datasets gives the 
tightest constraints, with wq = — 1.03±0.10, Wa ~ 0.30lo H 
and fim = 0.262 ±0.013. These values correspond to a lower 
pivot redshift Zp = 0.28 for which we get uide(zp = 0.28) = 
— 0.969±0.049. Even allowing for dynamic dark energy mod- 
els current observations can constrain the present value of 
iode at the 10% level. 

Combining the shape of the LRG correlation function 
with the rBAO measurements we obtain the constraints 
Wo ~ — l.OGlo lg and Wa = O.Ol? ?. These constraints give 
a much lower pivot redshift, with — 0.11, for which we 
find WBsizp = 0.11) = -1.06 ± 0.19. This shows that cur- 
rent large scale structure data allows for the present value 
of a redshift-dependent dark energy equation of state to be 
determined at the 20% level. 
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Figure 12. The marginalized posterior likelihood in the f2m — 
Cde plane for the ACDM parameter set extended with Hj. 
(Eq. l(9]|). The dashed lines show the 68 and 95 per cent con- 
tours obtained using CMB information alone. The solid contours 
correspond to the results obtained from the combination of CMB 
with the shape of the LRG The thick solid line corresponds 

to flat models, where Q/. = 0. The filled contours correspond to 
the 68% CL in each case. 

5.4 Non-flat models 

The flatness of the Universe is one of the generic predic- 
tions of most common models of inflation. A detection of a 
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Table 5. The marginalized 68% interval constraints on cosmological parameters for non-flat ACDM models (i.e. the parameter set 
defined by Eq. JSJl), obtained using different combinations of the datasets described in Section [2] as stated in the column headings. 





CMB 


CMB + i{s) 


CMB + rBAO 


CMB + SN 


PMR -1- 
+ rBAO 


PMB -1- ^f-j"! 
+ rBAO + SN 


i(s)+ rBAO 




-0.064+«'0f3 

— U.Uuo 


-0.0016+O-S 


-0.0036+«0053 


-0.012+0-010 


-0.0034+0-0051 

— U.UUOi 


-0.0035+0-0050 

— U.UUoi 


r, -19+0.089 
— U.Uoo 


100 e 


1 0412+°°°22 


1.0415+0-«E;^2 


1-0415+Hg22 


1 0414+0-0022 

— U .UU-iO 


l-0415+°-r22 


1-0415+0-0023 


i.4oiO:33i 




— U .UUDU 


0.1064+0-0«« 

— U .UU4o 


0.1086+0|;°1^ 

— U .UU'iO 


0.1080+0-°gS 

— U .UU4o 


0.1072l0:00« 


0.1084+0-0011 


o-ii2li::EJ^^ 


100 UJ^ 






2.272«;0f, 


2.26610-050 


9 97/1+0.050 

— 050 


2.267115:050 


3-itl:6 


T 


0.087lO;«l? 


0.089t°;EJ}? 


O.OQOlooiJ 


0-090l°:EJl^ 


0.089lO:Ul? 


0.08810:0^6 




ris 






0.963100}^ 


0-960l°:EJl^ 


0.963lO:«l^ 


0.96110:01^ 


0-99-0.31 




3.06it°;«|? 


o nc:9+0.038 

O.U0Z_Q ggy 


o.uoz_Q 039 


q nc:o+0.038 
o.UOO_o 038 


3 055+0-037 
O.UOO_o 038 


q nc:c:+0.035 
o.UOO_Q 035 


- 




0.56«:1? 


7P;f;+0.021 
U.COD_(, 


740+0 018 


n 717+0.023 

'J-' ^' -0.023 


0.748l°:«- 


741+0.013 
'^•'^1-0.013 


894+0-091 

U.O»'4_Q 090 




0.4981°:^ 


U.z^lo_Q 024 


U.ZDO_o 018 


205+0 030 
u.zyo_o 030 


0.256lO:«l« 


262+0-013 
u.zuz_o 013 


0.226lO:0« 




0.765i°:°g 


n 770+0. 028 
'J-' '«-0.028 


701+0.028 
'^•^ ^-^-0.028 


n 709+0. 028 
'J- '''^-0.028 


n 709+0.027 

'J-'*'^-0.027 


n 707+0.026 

O-'»'-0.026 




to/Gyr 


i5.7il:« 


-,0 71+0.39 
±0. ( i_o.38 


,0 oc: + 0.26 
10.00_o 26 


1 . 27+0-46 


13.8310:^5 


13.8510:25 


14.21^:0 




10.311:1 


10.41-3 


10.511:1 


10.4911:1 


10.411:^ 


10.311:^ 




h 




0.730t°:«i^ 


0.707lO:lil 


0.668lO:«|? 


r, 71 q+0.019 

U. / iJ_o.oi9 


707+0-017 
'J-"J'-0.017 


0.7810:15 



non-negligible curvature would have profound implications 
for our understanding of the mechanism thought to be re- 
sponsible for seeding density fluctuations in the Universe. 
The flatness hypothesis is also important because it has a 
strong impact on the constraints on the remaining cosmo- 
logical parameters, since many of them are degenerate with 
Q.k. In this section, we analyse non-flat models and include 
f2fc in our parameter space, assuming that the dark energy 
is given by vacuum energy with iude = — 1- 

When Q.k is allowed to float, the CMB data alone are 
unable to constrain all the parameters at the same time, 
giving rise to the well known geometrical degeneracy, which 
is completely analogous to the one between wde and Q.m 
described in Section[5]2] For each choice of Wb, Wdm and ilm, 
it is possible to find a value of fife (or JIde) which will give 
the same value of I a- This strong degeneracy can be seen in 
the dashed lines of Fig. 1121 which shows the two dimensional 
marginalized constraints in the plane J7ni — ^de- We plot 
the marginalized constraints on Sl^ in Fig. 1131 The dashed 
curve shows the results for CMB data alone. This narrow 
degeneracy produces poor marginalized constraints on the 
curvature of the Universe, with Q.k = — 0.064lo-063. 

Adding the independent constraint from the shape of 
the LRG ^(s) helps to break this degeneracy. The solid 
lines in Fig. [T^] show how with this extra information the 
two dimensional constraints in the plane fim — Ode close 
up over the locus of the flat models (shown by the thick 
solid line). From CMB data plus the LRG ^(s) we get 



fit 



-0.0016: 



Combining the CMB data with the 



other external datasets yields similar results, with Q,k 



-0 012+0-010 

U.UiZ_o 010 



-0.0036 



+0.0053 



for CMB plus SN and ^.k 
ft-om CMB plus rBAO. 

The combination of the shape of the LRG ^(s) and 
the rBAO data is not able to give meaningful constraints 
on fife. This data combination allows for a wide range of 
closed models giving a very poor one dimensional marginal- 
ized constraint of Q,k = —0. 12111:088 • On the other hand. 



CMB only 

CMB + f(s) 

CMB + f(s) + 
rBAO + SN 



0.4 
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Figure 13. The marginalized posterior likelihood in for the 
ACDM parameter set extended with Hfe (Eq. Q). The dashed line 
shows the likelihood obtained using CMB information alone. The 
solid line corresponds to the result obtained from the combination 
of CMB data with the shape of the LRG ^(s). The dot-dashed 
line shows the likelihood when all datasets are considered. 



these datasets can give a constraint on the matter density 
of fini = O.226I0 0I7, with a similar precision to the one 
obtained in the parameter spaces analysed in the previous 
sections. 

Our tightest constraint comes from the combination 
of the four datasets, which gives fife = — O.OOSSlo'oosi ^ 
which is dominated by the combination of CMB and rBAO 
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Figure 14. The marginalized posterior likelihood in the Qj, — 
WuE plane for the ACDM parameter set extended with and 
'"DB (Eq. Q). The dashed lines show the 68 and 95 per cent con- 
tours obtained using CMB information alone. The solid contours 
correspond to the results obtained from the combination of CMB 
data with the shape of the LRG ^{s). The dot-dashed lines show 
the constraints obtained using all data sets. The filled contours 
correspond to the 68% CL in each case. 

data. This resu l t is in full agreement with that found by 
iKomatsu et al. I l|2009l ) from the combination of WMAP-5yr 
data with the BAO measurements from Percival et al. of 
fife = -0.0050lo:oo60- This means that, with the hypothesis 
that vacuum energy is the source of the accelerated expan- 
sion of the Universe, current data allow us to probe spatial 
curvature up to fik — 5 x 10^'^. As we shall see in Sec- 
tion (53] this limit does not change when the assumption of 
wde = — 1 is relaxed. 



5.5 Dark energy and curvature 

The assumption of a flat Universe has important implica- 
tions for the constraints on the dark energy equation of 
state. Neither the curvature nor the dark energy play a role 
in determining the physical scales of the acoustic peaks in 
the CMB. These quantities do, however, have an impact on 
the angular diameter distance to the decoupling epoch. Fur- 
thermore, their contribution is degenerate. An increment in 
the value of flk, which reduces the value of Da (2*), can 
be compensated for by an increase in wde, which produces 
the opposite effect. This means that, instead of a one dimen- 
sional degeneracy as in the cases analysed in Sections l5.2l and 
15.41 this parameter space has an extra degree of freedom for 
every choice of Ub, <^dm and ilm- Fig. ll4l shows the marginal- 
ized posterior likelihood in the wde — plane. Due to this 
degeneracy, the region allowed by the CMB data, shown 
by the dashed lines, covers a wide range of this parameter 
space. When the information from the shape of the correla- 
tion function is included in the analysis, the extra distance 



measurement reduces this degeneracy to an (approximately) 
one dimensional region in the kjde — plane shown by the 
solid lines, but is also unable to place meaningful constraints 
on these parameters. 

An extra piece of information, coming from rBAO or SN 
data can break the remaining degeneracy. The dot-dashed 
lines in Fig. [14] correspond to the results obtained by com- 
bining the CMB measurements, the SDSS ^(s), rBAO and 
SN data. Combining all these datasets, the contours close-up 
over the canonical ACDM values, with Qk = -0.0018lo:oo54 
and Wde = — O.QBSlg Qgg It is important to note that the 
precision of these constraints is similar to those obtained 
with the assumptions of a flat Universe (Section 15. 2p or a 
cosmological constant (Section 15. 4|) . showing the robustness 
of our results. 

As in Sec. 15.41 the combination of the LRG ^(s) and the 
rBAO measurements is not able to give useful constraints in 
this parameter space. In this case, there is a wide range of 
closed models allowed by the data, which leads to a poor 
marginalized constraint of Qk ~ ~0-WQ^q'qq2- other 
hand, the large scale structure data is able to constrain the 
dark energy equation of state to wde ~ —0.91 ± 0.14, with 
the same level of accuracy obtained with the assumption of 
a flat Universe. 



6 DISTANCE CONSTRAINTS 

In this section we explore the distance measurements that 
can be obtained from the shape of the LRG redshift space 
correlation function and its combination with other datasets. 
For the results of this section we have imposed the conser- 
vati ve constraint on t h e LRG bias factor of 6 < 4. 

lEisenstein et al.l (|2005h used the correlation function 
of a sample of around 46,000 LRGs from the SDSS to ob- 
tain constraints on fl^nh^ and the combinations of distance 
measurements Dv{z) (defined in Eq. [16} and A{z), defined 
as 

A(z)^Dv{z)^^^^, (26) 
zc 

at the mean redshift of the survey Zm ~ 0.35 and found 
-Dv(zm = 0.35) = 1370 ± 6 4Mpc and A(zm = 0.3 5) = 
0.469±0.017. In their analysis. lEisenstein et aTl (120051 ) fixed 
the values of Wb ~ 0.0223 and ris = 0.98, as well as adopt- 
ing a fixed BAO damping scale fe* = O.llftMpc"^ and 
varied only Qmh^ and i3v(2m = 0.35). Fixing the values 
of these parameters has implications for the obtained con- 
straints. They also implemented a different model for the 
non-linear distortion of the shape of the correlation func- 
tion and included a wider range of scales in their anal- 
ysis (10/t~^Mpc < r < 177/i~^Mpc) than we consider. 
lEisenstein et al.l l|2005h used a fixed template for the nonlin- 
ear distortion of the linear perturbation theory correlation 
function, ignoring changes expected on varying the cosmo- 
logical parameters (e.g. as). Also, no attempt was made to 
model scale dependent effects such as redshift space distor- 
tions or bias, which could become important for pair separa- 
tions approaching 10h~^ Mpc. In our analysis, we only con- 
sider pair separations given by r > 42.5 /!-~^Mpc, for which 
our theoretical model is in excellent agreement with the re- 
sults of N-body simulations (see Fig. [3]). For the wide priors 
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Table 6. The marginalized 68% interval constraints on the cosmological parameters allowing for non-flat models and variations in 
the (redshift independent) dark energy equation of state parameter (i.e. the parameter space defined by Eq. JSj), obtained using 
different combinations of the datasets described in Section [21 as stated in the column headings. 





CMB 


CMB 4- F(<i) 


CMB + rBAO 


CMB + SN 


CMB + 
+ rBAO 


CMB + 
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Figure 15. The marginalized posterior likelihood in the 
DAizm = 0.35) - H(2in = 0.35) plane for the ACDM parame- 
ter set extended with Oj. and wde (Eq. lO). The short-dashed 
and long-dashed lines show the 68 and 95 per cent contours ob- 
tained individually from CMB information and the LRG ^(s) re- 
spectively. The solid contours correspond to the results obtained 
from the combination of CMB data plus the shape of the LRG 
5(s), which follow approximately a line of constant G(2m = 0.35) 
(dotted line) defined by Eq. II28I I. The dot-dashed lines show the 
contours obtained including also rBAO information. The shading 
fills the 68% CL contours. 



listed in Table 1, and using only information from the shape 
of ^(s) over the range of pair separations 42.5 h~^Mpc < r < 
150 /i^^Mpc, we find weaker constraints, Dv(-2m = 0.35) = 
1230±220 Mpc and A{z^ = 0.35) = 0.424±0.064. To achieve 
an accuracy comparable to that of Eisenstein et al. it is nec- 
essary to combine the ^(s) constraints with CMB informa- 
tion. In this case we find Dv(-Zm = 0.35) — 1300 ± 31 Mpc 
and A{zrn = 0.35) = 0.447 ± 0.015. 

In the analysis o f the fi fth year data of the WMAP satel- 
lite , iKomatsu et al. ] l|2009l ) produced a set of distance priors 
that contain most of the information in the WMAP power 
spectrum. This set contains (i) the physical baryon density 
w^j, [a) the redshift to the decoupling epoch z,, (tii) the 
"acoustic scale" of Eq. p9|) and {iv) the "shift parameter" 
R given by 



R{z,) = DAiz.) 



(27) 



with their respective covariance matrix. Here, we provide 
an extended set of distance priors combining the information 
from the clustering of the LRGs with CMB data. For this, we 
expan d the set of parameters provided by iKomatsu et al. I 
(120091 ) by adding an extra constraint on the quantity 



G{Z^) = DAiZm) X h{z^y 



(28) 



where h{zin) ~ i/(2:m)/100kms~* Mpc~^. Now we show 
that this parameter is well constrained by the combination 
of the CMB measurements and the LRG ^(-8). 

Fig. [15] shows the marginalized posterior likelihood in 
the DA(zn-i) — H{zm) plane for the full parameter space of 
Eq. ((3)1. The constraints obtained individually from CMB 
information (short-dashed lines) and the LRG (,{s) (long- 
dashed lines) show strong degeneracies between these pa- 
rameters. The combination of the two datasets reduces the 
allowed region for these parameters to a narrow degeneracy 
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Figure 16. The marginalized posterior likelihood in the iode ^ 
Qui plane for the ACDM parameter set extended with iode 
(Eq. 10). The dashed lines show the 68 and 95 per cent contours 
obtained using the full combination of CMB data plus the shape 
of the LRG The solid contours correspond to the results ob- 
tained using the set of extended distance priors of Tables Bl and 
B2 of Appendix B. The shading fills the 68% CL contours. 



(solid lines) which follows approximately a line of constant 
G (shown by the dotted line). For this combination we find 
the constraint G{zm = 0.35) = UTSti? Mpc (a 1-4% error) 
which we include in our set of extended distance priors. The 
degeneracy between Da and H can be broken by including 
an extra piece of information. In this case it is possible to 
obtain separate constraints on these parameters. Including 
the rBAO measurements (as shown by the dot-dashed lines 
in Fig. [151) we obtain DA{zm = 0.35) = 1363 ± 34 Mpc and 
H{zru = 0.35) = 83.3 ± 2.2 km s"^ Mpc'^ By using the SN 
data instead we obtain DA(zm = 0.35) = 13361^!^ Mpc and 
H{zm = 0.35) = 86.4±5.6kms-iMpc-\ 

The mean values and covariance matrix of the set of 
extended distance priors obtained from the combination of 
CMB information and the shape of the LRG ^(s) are listed 
in Tables Bl and B2 of Appendix B. Fig. ll6l shows a compar- 
ison of the marginalized constraints in the plane flm — wde 
for the parameter space of Eq. (O obtained by using the 
set of extended distance priors (solid lines) and the ones ob- 
tained using the full CMB and ^(s) data (dashed lines). This 
shows that this set of constraints contains the most relevant 
information from the combination of these datasets and can 
be used to replace them to obtain constraints on cosmo- 
logical parameters in combination with other cosmological 
observations like SN data or weak lensing, simplifying and 
accelerating the required numerical procedure. 



7 CONCLUSIONS 

In this paper, we have applied a new model for the shape 
of the correlation function at large pair separations to 



th e measurement of LRG clust ering from the SDSS made 
bv ICabre fc Gaztafiaga I (|2009al ). Using a fuU Monte-Carlo 
markov chain analysis and combining the LRG measure- 
ments with the latest compilations of CMB and SNe data, 
we have presented a comprehensive set of constraints on cos- 
mological parameters for different combinations of datasets 
and for different parameter spaces. 

Large scale structure measurements and modelling have 
now reached a level of precision where they provide con- 
straints on cosmological parameters which are competitive 
with those obtained from other datasets. Using only the 
LRG correlation function and the measurement of the radial 
BAO peak, it is possible to determine the dark energy equa- 



tion of state parameter to uide = — 1.05_ 



This provides 



further confirmation of the dark energy scenario, indepen- 
dent of CMB or SN observations. 

The availability of several high quality datasets gives 
us the opportunity to scrutinize the derived constraints and 
to isolate possible systematic effects in either the observa- 
tions themselves or the theoretical model used to describe 
them. In general, different datasets are sensitive to different 
parameter combinations and the constraints are optimized 
when datasets are combined. However, in the case in which 
datasets are responsive to similar parameter combinations, 
it is important to check that consistent results are obtained 
before combining the datasets. It would be meaningless to 
combine two such datasets if a tension existed between the 
parameters values returned from each one individually. 

The tables of parameter constraints given in this pa- 
per show that, to first order at least, the "Union" super- 
nova dataset of iKowalski et al. I (|2008l ) and the LRG clus- 
tering from the SDSS yield very similar constraints on cos- 
mological parameters when analysed in combination with 
CMB data, for all of the parameter spaces considered. As 
shown in Section 15.21 the recovered value of the dark en- 
ergy equation of state in the CMB plus ^(s) case of u;de ~ 
—0.996 ±0.090 shows a remarkable agreement with that ob- 
tained from the combination of CMB with the SN data of 
uiDE ~ —0.950 ± 0.055. The consistency of the constraints 
we obtain from different combinations of datasets is a reas- 
suri ng validation of our an alysis technique. 

iPercival et al. I (l2007d ). however, found a tension at the 
2.4(7 level between the constraints on Dv{z = 0.35)/ Dv{z — 
0.20) coming from the BAO signal measured from galaxy 
samples drawn from the 2dFGRS and the SDSS (includ- 
ing a smaller LRG sample t han the one considered here), 
and the SN data from Asti er et al. I (2006). When folded 
with the CMB data, this discrepancy leads to different pre- 
ferred values for the dark energy equation of state, with the 
SN data preferring uide — —1 and the BAO data point- 
ing to wuE < —1 with a significa nce o f 1.4ff (see their fig- 
ures 12 and 13; see also Lazkoz et al. I l|2008l )). Our analy- 
sis shows that the LRG clustering and SNe data give con- 
sistent constraints on the distance scale Dy jz = 0.35). 
The tension reported bv lPercival et al. I (|2007ch could result 
from the BAO signal at z = 0.2 (which corresponds to the 
main 2dFGRS and SDSS galaxy samples). This is unlikely 
to be the answer, however, sin ce the rBAO measurements 
from iGaztafiaga et al71 l|2008bl ) , which separate the signal 
at z = 0.24 and z = 0.43, also give a consistent answer with 
uiDE ~ —0.92 ±0.15. The treatmen t of the BAO feature is 
quite different in the two analyses. iPercival et al~l (|2007d ) 
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effectively take a BAO template generated for a given set 
of cosmological parameters and then fit this to the mea- 
sured BAO features for different parameters. Furthermore, 
the technique of Percival e t al. requires a referenc e spectrum 
to be defined. Finally, as iPercival et al. I l|2007cl ) speculate 
themselves, the problem could lie in the way in which the 
BAO are damped; these authors assume a fixed damping 
scale, whereas we treat the damping scale as a parameter. 

One criticism that could be levelled at our analysis 
is the apparent level of agreement between our new theo- 
retical model and the measured LRG clustering. We have 
shown that the model is remarkably accurate when com- 
pared to the correlation function measured for the dark 
matter or for s amples of haloes drawn from N-body simula- 
tions (see also lCrocce fc Scoccimarro 11200 ISanchez et al. 1 
I2OO8I ). However, in the case of the observed LRG clus- 
tering, the relative height of the BAO feature is greater 
than predicted by the model. This is illustrated in Fig. 1171 
This behaviour is not pecu liar to the measurement by 
ICabre fc Gaztanaga I (l2009ah. The e a rlier m easurement of 
LRG clusteri ng by Eisenstein et al. I |200^ had a similar 
form (see also lOkuniura et al. II2OO8I ). A key result from this 
paper is the dependence of the parameter constraints on the 
minimum pair separation used in the analysis (Fig. 5). The 
mean parameter value returned is insensitive to the choice of 
the minimum pair separation; the only noticeable difference 
is an increase in the Icr range on selected parameters as fewer 
data points are included. Similar conclusi ons were reached 
using a less challenging range of scales by lOkumura et al. I 
l|2008l l. This confirms that the information we extract from 
the overall shape of the correlation function is consistent 
with that contained in the BAO feature. 

Fig. [17] also considers an extension to our best fitting 
model which involves allowing for a constant shift in ampli- 
tude (Sec 2.1.1). We have tested that the marginalization 
over a constant additive term has a negligible impact on the 
cosmological constraints that we derive, which are driven by 
relative changes in the goodness of fit. However, the addi- 
tion of this extra nuisance parameter can improve the ap- 
pearance of the fit and therefore reduce the absolute value of 
■ If we use an absolute fit with the full covariance ma- 
trix, we find a probability of about 25% that our best fitting 
model describes the measured correlation function (the ex- 
act value varies between 20% and 33% depending on which 
parameters we consider as degrees of freedom). If we allow 
for a constant systematic shift, as argued in section [2. 1.1 1 the 
best fit model for the cosmological parameters remains the 
same, but the probability of the model describing the data 
improves slightly to 60% in the range K = 0.003 — 0.006. 
This improvement is not significant enough to justify the 
extra parameter, despite the visual impact of this shift in 
Fig. El 

We used the combined information of the CMB mea- 
surements and the shape of the LRG f (s) to obtain a set of 
distance priors that combine the most relevant information 
from these data sets to place constraints on dark energy 
models and spatial curvature. For this, we expanded the 
set of distance priors provided by iKomatsu et aTl (|2009l ) 
by adding a constraint on the quantity G{zni), defined in 
Eq. (|28|l . which is well constrained by the combination of 
these datasets (with an error of less than 2%). The mean 
values and covariance matrix of this extended set of dis- 
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Figure 17. Circles with error bars show t he redshift space cor- 
rclation function used in this paper (from ICabre &: Gaztaiiaga I 
The solid line corresponds to our best-fit model within 
the framework of the basic ACDM parameter space described in 
Sec. 15.11 while the dashed line includes in addition a systematic 
K-shift with K = 0.0045. This shift may originate from unknown 
systematic effects as described in Sec. 12.1.11 Even though the vi- 
sual impact seems important, this constant shift (note the bias 
factor also changes) does not compromise the goodness of fit of 
our best fitting model or the range of cosmological parameters 
derived in this paper. Our best fitting model (ie using K = 0) 
has a probability of 24 % to be in agreement with the data, as 
compared with 62 % for the model with K = 0.0045. Thus there 
is not a significant improvement on adding a constant shift to the 
model. 



tance priors are listed in Tables Bl and B2 of Appendix 
B. These can be used instead of the combination of CMB 
and 5(s) data used in our analysis to obtain constraints on 
cosmological parameters in combination with other cosmo- 
logical observations, making the required calculations faster 
and simpler. 

Our analysis shows that current large-scale structure 
datasets (using the shape of the correlation function at 
large pair separations, along with the form of the BAO 
peak) give constraints on the basic cosmological parame- 
ters that are competitive with those from the best available 
SNe datasets. Moreover, possible systematic distortions in 
the clustering signal have been modelled extensively. The 
theoretical model applied in this paper has been shown to 
provide an accurate description of the clustering for the level 
of error expected in volumes up to two orders of magnitude 
larger than the current LRG datasets. This bodes well for 
experiments planned for ten year's time, which will produce 
surveys covering vol umes in the region of 100/i-^Gpc^ (e.g. 
ICimatti et al. I2OO8I ). Nevertheless, it will be prudent to pur- 
sue as many approaches as possible to measure dark energy 
to allow comparisons to be made between the results ob- 
tained with different datasets, as carried out in this paper. 
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APPENDIX A: THE THEORETICAL 
MOTIVATION FOR THE MODEL FOR THE 
SHAPE OF i{R) 

In recent years substantial progress has been made in 



density fluctuations ( 


Crocce & Scoccimarro 1 l2006al. l2008l: 


Matarrese & Pietroni 


I2OO7I: iMatsubara 20081: iPietroni 1 


2008'; Taruva & Hiramatsu 112008^. Within the theoretical 



framework of Renormalized P erturbation Theory (RPT), 
ICrocce fc Scoccimarro I l|2006al lb[) showed that the evolution 
of the power spectrum from its initial value in the linear 
regime, Ph{k), can be described as the sum of two contribu- 
tions 



Pnl (fc, z) = Pl (fc)G2 {k,z)+PMc{k,z), 



(Al) 



where G is the nonlinear propagator which weights how 
much power can be directly linked to the linear epoch (e.g. 
damping the higher BAO harmonics), and Pmc represents 
the new variance generated by mode coupling. 

In configuration space, the first term in Eq. (|A1|I leads 
to a convolution of the linear theory correlation func- 
tion with a nearly Gaussian k ernel (since G ~ e"*" , see 
ICrocce fc Scoccimarro I (|2006lJ l). which cause s the acoustic 
peak to broaden and shift to smaller values l|Smith et al. I 
2008; Cr occe fc Scoc cimarro 2008). In turn, the contri- 
bution of the mode-coupling term can be computed to 
leading order at large pair separation, yielding ^Mc{r) oc 
^i{r)^i^\r) (|Crocce fc Scoccimarro I I2OO8I ). where ^i{r) is 
the first derivative of the linear theory correlation function 
and ^(r) is defined by Eq ([TT}. 

In fact, Eg.f IAl| l is also valid for the nonlinear spectra 
of density (5)/velocity divergence (6^) fields, Pge and Poe, 



PNL,ab = PhGaGb + PMC,ab, 



(A2) 



with a — 5,6. 

To study the impact of redshift distortions and bias we 
put together the di scussion abo ve and the ansatz given in 
IScoccimarrol (120041 ') (see also Percival fc Wh ite (2009)); 



P.(k) = e-'='^'^'"" [b^Pss + 2ft^bPse + f^i''Pe, 



(A3) 



where Pss , Pse and Pgg were introduced above and a local 
linear bias relation between the density fluctuations in the 
distribution of galaxies and dark matter is assumed (but 
with unbiased galaxy velocities). The damping factor in 
Eq. ()A3|) arises after assuming a Gaussian PDF for the pair- 
wise velocities with rms (if evaluated in the linear regime). 
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- 4 / 

In Eq. (|A3P / is the logarithmic derivative of the linear 
growth rate with respect to the scale factor, f — d InD/d Ina, 
and n = kz/k is the cosine of the line-of-sight angle. Note 
that the small scale redshift space distortions, commonly re- 
ferred to as "fingers of God" (FOG) due to the elongation of 
virialised structures, are difficult to model accurately, partly 
due to the non- Gaussianity of the PDF, which persists even 
to large scales (jScoccimarro Ifiool ). In addition, the damp- 
ing factor depends on galaxy type (or mean halo mass and 
satellite population) 

We are interested in the monopole of the correlation 
function, i.e. the angle average of the Fourier transform of 
Pa{^) in Eq. (|A3|) . Since the time required to evaluate the 
different power spectra is not well suited to the exploration 
of large parameter spaces, we would like to motivate a para- 
met ric description of the pr o blem. 

ICrocce fc Scoccimarro I (|2006lJ ) showed that both Ge 
and Gs are of approximately Gaussian form, and with sim- 
ilar smoothing length. Therefore the main contribution to 
the correlation function at large separations can be mod- 
eled, combining Eqs. (|A3IA2p . as oc ® e"^''*'"'^ The FOG 
prefactor in Eq. (|A3|l leads to an extra suppression, which 
is subdominant to that of nonlinear gravity for low redshift 
data (where / ~ 0.5). However, it can be accounted for to 
some degree by smaller values of fc^. 

The mode c oupling spectra in E q. (|A2[| do differ at the 
several % level (|Scoccimarro ][200i), but they aU give con- 
tributions scaling as ~ (the derivative of the correlation 
function) in configuration space at large separations. The 
convolution with the FOG suppression leaves this contribu- 
tion unaltered as the smoothing length introduced by the 
FOG is smaller than features present in even at BAO 
scales. 

Therefore we are left with the following parametriza- 

tion, 

^NL(r) = b^{Ur) ® e-<'=*'-)' + ^mc ?L(r) ^i'\r)}, (A5) 

for the correlation function at large pair separations in red- 
shift space, where b and ^mc are nuisance parameters to 
account for bias and also the enhanced correlation ampli- 
tude due to redshift distortions (Kaiser 1987). 

The assumptions behind the model in Eq. (|A5|I . one 
could argue, are over-simplifications, particularly for biased 
tracers such as LRGs (Zehavi et al. 2005; Blake et al. 2007), 
or for deeper surveys where a more systematic derivation of 
redshift distortions may be needed. Nonetheless this form 
is well motivated given the size of statistical errors in the 
SDSS-D R6 LRG sample , and based on results from Sec. 13.21 
(see also lSanchez et al. I (|2008l )) should be accurate also for 
future larger surveys. 



Table Bl. The mean values and variances of the parameters 
in the set of extended distance priors described in Section |6l us- 
ing the compilation of CMB data in combination with the LRG 
correlation function. 



Parameter value 



100<.^b 2.28 ± 0.55 

2* 1090.12 ±0.93 

£a(2*) 301.58 ±0.67 

R(z*) 1.701 ±0.018 
G{zm) 1175 ±21 



can be used to replace the CMB and ^(s) data in multipara- 
metric analysis since it contains most of the information of 
the combination of these two datasets. With this simplifi- 
cation, the likelihood of a given cosmological model can be 
easily computed as -2 In £ oc (D - T)' C"^ (D - T)* where 
D is a vector containing the constraints of Table Bl, T con- 
tains the corresponding prediction for a given cosmological 
model and is the inverse of the covariance matrix of Ta- 
ble B2. This procedure considerably reduces the computing 
time required to obtain constraints on cosmological param- 
eters compared with the full analysis of the data and can be 
used in combination with other datasets like weak lensing 
or the Hubble diagram of type la supernovae. 



APPENDIX B: EXTENDED DISTANCE 
PRIORS 

Here we list the constraints on the set of extended distance 
priors described in Section [S] Table Bl lists the mean val- 
ues and variances of the parameters in this set and Table B2 
gives the corresponding covariance matrix. This information 
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Table B2. The covariance matrix of the parameters in the set of extended distance priors described in Section |S] 



lOOi^b z, ^a(2.) Riz*) G{z^) 

100a;b 2.99 X 10^^ -4.2030 X 10"* -1.9988 X 10"* -3.4393 x -1.9978 x IQ-^ 

z, 8.6812 X IQ-i 0.2557 1.1999 x IQ-^ 3.2355 

£a(^.) 0.4558 3.1265 x 10"^ 1.8247 

R{zt) 3.1460 X 10"* 3.9649 x 10"^ 
G(2m) 4.3784 X 10^ 
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